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Abstract: The holographic dual of large- A^^. super- Yang- Mills coupled to a small number 
of flavours of fundamental matter, N( <^ N^, is described by probe D7-branes in the 
gravitational background of A^^ black D3-branes. This system undergoes a first order phase 
transition characterised by the 'melting' of the mesons. We study the high temperature phase 
in which the D7-branes extend through the black hole horizon. In this phase, we compute 
the spectral function for vector, scalar and pseudoscalar modes on the D7-brane probe. We 
also compute the diffusion constant for the flavour currents. 
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1. Introduction 



The gauge/gravity duality has created a powerful framework for the study of strongly coupled 
gauge theories [1, 2, 3, 4]. Thermal properties of such theories are of considerable interest both 
in their own right and in connection with the experimental program on heavy ion collisions 
at RHIC and LHC. 

In this paper, we study certain aspects of finite-temperature behaviour of strongly coupled 
N = 2 super- Yang-Mills theory with dynamical quarks in Minkowski space. The holographic 
description of the theory with gauge group U{N^) and hypermultiplets in the fundamental 
representation (as well as one adjoint hypermultiplet) is given by the well studied system of 
D3- and D7-branes. In the limit of large and large 't Hooft coupling A = gy^^c with fixed 
A^f, the theory is described by Ni probe D7-branes in the near-horizon geometry of D3- 
branes, i.e., AdSs x [5]. At finite temperature, the background geometry contains a black 
hole [6] . Although the M =1 theory is non-confining and thus no confinement-deconfinement 
phase transition is expected at finite temperature, the presence of the quark mass scale Mq 
leads to a first order phase transition for fundamental matter [7]. The transition occurs at a 
temperature T ~ Mq/\/\ and is characterized by the dissociation or 'melting' of the mesons, 
i.e., the quark-antiquark bound states^. 

Our goal is to study thermal dissociation of the bound states as well as the flavour 
current diffusion by computing spectral functions of mesonic operators in the framework of 
gauge/gravity duality. 



Minkowski embedding Criticai embedding Black hoie embedding 




Figure 1: Various possible D7-brane embeddings in the black D3-brane geometry with temperature 
increasing from left to right. 

The thermal behaviour of the M =1 gauge theory with fundamental matter has recently 
been actively studied [8, 9, 7, 10]. The holographic description gives a simple picture of the 
phase transition. If the D7-branes are sufficiently far from the event horizon of the black hole, 
they are gravitationally attracted towards the horizon but their tension is sufficient to balance 
this attractive force. The probe branes then lie entirely outside of the black hole in what 
were denoted 'Minkowski' embeddings in [7, 10] - see fig. 1. As the temperature is raised, 
both the radial position and the energy density of the event horizon increase. Therefore, 

^Recall that for these supersymmetric field theories, the fundamental matter includes both fermions and 
scalars, which we will refer to collectively as 'quarks'. 
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above some critical temperature, the gravitational attraction of the black hole overcomes the 
brane tension and the D7-branes are pulled into the horizon. These configurations where 
the branes fall through the horizon are referred to as 'black hole' embeddings. In between 
these two branches, there also exists a critical solution which just 'touches' the horizon. 
Thermodynamic considerations show the latter is by-passed by a first order phase transition, 
in which the probe branes jump discontinuously from a Minkowski to a black hole embedding. 

In the dual field theory, this phase transition is exemplified by discontinuities in various 
physical quantities, e.g., the quark condensate. However, the most striking feature of the 
transition is found in the spectrum of the mesons. The latter correspond to excitations sup- 
ported on the probe branes - see, e.g., [11, 12]. In the low-temperature or Minkowski phase, 
the mesons are stable (to leading order within the approximations above) and the spectrum 
is discrete with a finite mass gap. In the high-temperature or black hole phase, mesons are all 
destabilized and rather one finds a continuous and gapless spectrum of excitations^ . Accord- 
ingly, spectral functions of mesonic operators in the low-temperature phase are characterized 
by (5-function peaks^ (with the decay width of these particles and the continuum contribution 
of multi-particle states both suppressed by factors of 1/N^), whereas in the high-temperature 
phase spectral functions are essentially featureless^ (fig. 2). More interesting behaviour is 
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Figure 2: Sketch of typical spectral functions in the (a) low-temperature (Minkowski) and (b) high- 
temperature (black hole) phases. 



observed when the system evolves from the high-temperature phase into the low-temperature 
phase through the metastable 'supercooled' phase. We show that in this case the serene land- 
scape of fig. 2b is distorted by peaks corresponding to quasiparticle excitations, and these 
excitations are eventually transformed into the stable resonances shown in fig. 2a. 

These features of the spectral functions are controlled by the analytic structure of the 
corresponding retarded correlators in the complex frequency plane. In the high-temperature 
phase, the poles of the retarded correlators (with the exception of the poles corresponding to 



^For a discussion of thermal resonances in the context of holographic models see e.g., [13]. 

derivation of the scalar meson spectral function at T = appears in appendix B. 
''The temperature-dependent part of the spectral functions exhibits damped oscillations with the period 
proportional to a Matsubara frequency, see section 2.3 for details. 
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hydrodynamic excitations) are located at a finite distance from the real axis thus making the 
spectral function featureless. As the temperature is lowered with respect to the quark mass 
scale, the poles move closer to the real axis and the spectral functions exhibit distinct peaks. 
Holographically, the poles of the retarded correlators correspond to quasinormal modes of the 
gravitational background [14, 15]. Numerical investigation of the full quasinormal spectrum 
of the D3-D7 system faces certain technical difficulties. In this paper we focus on computing 
the spectral functions for which the numerical methods are reliable. 

In addition to characterizing quasiparticle excitations of a thermal system, spectral func- 
tions also carry information about the medium's transport properties. Adapting techniques 
from [16, 17, 18], we compute the quark diffusion constant as a function of the parameter 
Mq/T in the high-temperature phase, and attempt to give a qualitative description of its 
dependence on the coupling for the full range of temperatures. 

Thermal dissociation of mesons as well as transport properties of the quark-gluon plasma 
can be studied in lattice QCD with the help of indirect methods such as the maximal entropy 
method [19, 20, 21, 22, 23, 24, 25, 26, 27, 28]. These studies suggest, in particular, that mesons 
survive as relatively well-defined resonances at temperatures well above Tc (2 — 3Tc). While 
the uncertainties of these lattice methods remain large, the holographic approach used in this 
paper serves as a source of quantitative and often analytically exact results for qualitatively 
similar finite-temperature models. 

An overview of the paper is as follows: in section 2 we review properties of thermal 
spectral functions in field theory and outline methods of computing spectral functions from 
dual gravity. These methods are illustrated by a simple example of computing the spectral 
function and diffusion constant for R currents in = 4 SYM. For vanishing spatial mo- 
mentum, the correlator, quasinormal spectrum, and the spectral function can be computed 
analytically. In section 3 we introduce the D3/D7-brane framework and review the D7-brane 
embeddings and thermodynamics. In section 4, we turn to the calculation of the spectral 
function for various mesonic operators in the high temperature phase of the M = 2 gauge 
theory. We consider a vector operator in section 4.1. In the special case of vanishing quark 
mass, we determine the spectral function analytically. In general, for arbitrary quark mass, 
the vector spectral function is computed numerically. In section 4.2, we turn to spectral 
functions for scalar and pseudoscalar operators, which are again bilinear in the fundamen- 
tal fields. Section 5 presents three independent computations of the diffusion constant for 
'light' quarks, using the membrane paradigm method [17], the Green-Kubo formula, and by 
calculating the lowest quasinormal frequency for the vector field on the D7-brane. Section 6 
contains discussion and observations about our results. Some details of our analysis are rele- 
gated to appendices: appendix A provides a review of the holographic dictionary relating the 
D7-brane worldvolume fields to their dual operators in the gauge theory; appendix B contains 
a derivation of the scalar spectral function at T = 0; appendix C provides a derivation of 
the high frequency asymptotics for the spectral functions; appendix D gives a partial analysis 
of the quasinormal modes for the pseudoscalar and scalar excitations; and finally, appendix 
E extends the computation of the quark diffusion constant, described in section 5.1, to the 
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holographic framework described by a Dq-brane probe in a near-extremal Dp-brane throat. 
2. Prelude: spectral functions and holography 

In general, finite-temperature correlation functions of conserved charge densities carry infor- 
mation about a medium's transport properties and quasiparticle excitations. This information 
is given, roughly, by the poles and the corresponding residues of the correlators, or, equiva- 
lently, by their spectral functions. Recently the study of these objects has been used to great 
effect in a holographic framework to study the thermal properties of various strongly coupled 
field theories [40]. In a holographic setting, the spectral functions are often easier to compute 
than the full correlators on the gravity side. According to the holographic dictionary, the 
poles are determined by the quasinormal spectrum of a dual bulk field fluctuation, whereas 
the spectral function is given by the imaginary part of the retarded correlator which is inde- 
pendent of the radial coordinate [15, 29]. In this section, we combine the necessary tools for 
computing the spectral functions from dual gravity and analyzing their properties, and then 
illustrate this technique using the simple example of strongly coupled A/" = 4 supersymmetric 
SU{N^) Yang-Mills (SYM) theory at large N^. In this case, the i?-current spectral function 
has been analyzed elsewhere [29, 30, 13] but we present a new analytic result (for vanishing 
spatial momentum). 

2.1 Field theory picture 

A thermal spectral function of an operator O is defined as^ 

^(a;,q) = J d'^x e--*+^'i-([0(i, x), 0(0)]) , (2.1) 

where the correlator is computed in thermal equilibrium at a temperature T. The spectral 
function is proportional to the imaginary part of the retarded correlator, 

^l(a;,q) = -2ImG^(a;,q), (2.2) 

where 

G^(a;,q) = -iJ d^x e-^'^^+^i^ e(x°)([0(t, x), 0(0)]) . (2.3) 

If O is an operator of a density of a conserved charge in a rotation invariant theory, the 
retarded thermal two-point function is determined by two independent scalar functions. In 
Fourier space, the correlator can be decomposed into the transverse and longitudinal parts 
[29] 

G^,{u, q) = Pj, n^(a;, q) + Pj^, n^(a;, q) , (2.4) 

''Our metric convention is (— , -h, +■, +)■ We assume translation invariance to be an unbroken symmetry of 
the theory. 
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where the index structure is absorbed into two mutually orthogonal projectors and P^^. 
Without loss of generality we can take the spatial momentum oriented along the direction, 
so that kn = {—u, 0, 0, q), with P = —uJ^ + . Then one has [29] 

G%A^) = G%A^)=^^^^^<1) ■ (2-5) 
Other components of the current-current correlator are 

G«W = -7237^n^(^'9)' (2-6) 
Glm = -^^^^\^.^). (2.7) 

G%A^) = -^^^^\^.^)- (2-8) 

In the long-time, long- wavelength limit (i.e., for ujjT ^ 1, qjT <^ 1) the functions ]l^{uj,q) 
and Il^{oj,q) have a universal behaviour dictated by hydrodynamics: IL^{(jj,q) is nonsingular 
as a function of the frequency, while Il^{io, q) has a simple pole at 

uj = -iDq^ , (2.9) 

where D is the charge diffusion constant. 

The rotation invariance implies that in the limit of vanishing spatial momentum at fixed 
a; > the two scalar functions coincide: 11-^ (a;, 0) = n^(a;, 0) = n(a;). Correspondingly, at 
q = one can define 

^(w) = •H^i^i(a;,0) = •H^2^2(a;,0) = 'H^^a^a (w, 0) . (2.10) 

The Green-Kubo formula relates the diffusion constant to the zero-frequency limit of the 
spectral function 'H(a;): 

limi-^(a;). (2.11) 

Here H is the charge susceptibility. The susceptibility is determined by the thermodynamics 
of the system in a grand canonical ensemble. 



„ (9n(/i) 
djjL 



(2.12) 

/i=0 



where n{^) is the charge density, /i is the corresponding chemical potential. 

In addition to hydrodynamic poles, the retarded correlators may have other singularities 
located in the lower half-plane of complex frequency. Assuming one of these singularities is a 
simple pole, 

qR- 1 



UJ — r2(g, a) -|- ir(g, a) ' 
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where a represents a set of parameters relevant for a particular theory, for the spectral function 
one has ^ 

Thus in the vicinity of w = 0, the spectral function has a peak characterized by a width ~ F 
and a height ('lifetime') ~ l/F. The peak has a quasiparticle interpretation if F ^ Q. 

The spectral function '^{lo) also has a characteristic form in the high frequency limit. 
This behaviour is determined by the leading short-distance singularity 

A 



lim J0(t,x)0(O)) = -.2|A +---^ (2-13) 



where A denotes the dimension of the operator O and ,4 is a dimensionless constant. A 
Fourier transform then leads to the following contribution to the spectral function 

^{u) ~ Au?^'"^ . (2.14) 

2.2 Gravity picture 

In the dual gravity picture, the conserved current couples to a boundary value of the gauge 
field fluctuation A^^ propagating in a specific gravitational background. One can form two 
gauge-invariant combinations of the fluctuation whose equations of motion (supplemented 
with appropriate boundary conditions) contain (in the limit where the gravity description is 
valid) full information about the functions n^(a;, g) and ^{uj^q) introduced in section 2.1. 
These gauge invariant combinations arc the transverse and longitudinal (with respect to a 
chosen direction of the spatial momentum) components El of the electric field in curved 
space [29] . Quasinormal spectra of the fluctuations Et and El determine the position of the 
poles of Ii^{uj, q) and n^(a;, q) in the complex u plane. 

2.3 A simple example: spectral function of R currents in = 4 SYM 

Correlators of /^-currents in strongly coupled M = 4 SU (N^) supersymetric Yang-Mills (SYM) 
theory at large were previously studied by means of the AdS/CFT correspondence both 
at zero [31, 32] and finite temperature [16, 33, 29, 30, 13]. 

In thermal J\f = 4 SYM, the retarded two-point correlators of the SU R-symmetry 
currents are determined by two independent scalar functions^, Il^{u!,q) and 11^ (a;, g). 

The holographic dual of thermal A^ = 4 SYM in fiat space is well known. The supergravity 
background describing the decoupling limit of black D3-branes is (see, e.g., [4]) 

ds' = ^ {-fir)dt' + c/x^) + ^ + -'dnl^ ' ^0123 = , (2.15) 

where /(r) = 1 — rg/r^ and the dilaton is constant. The horizon lies at r = tq and the radius 
of curvature L is defined in terms of the string coupling constant and the string length scale 



®In an equilibrium state without chemical potentials for the i?-charges, the correlation function of R currents 
has the form C^^ = d^^C^vix). In all expressions, the factor (5"'' is omitted. 
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4 as L'^ = 4:7rg^N^£f. According to the duality, originally proposed by Maldacena [1], type 
IIB string theory on these backgrounds is dual to four-dimensional J\f = 4 super- Yang-Mills 
(SYM) SU{N^) gauge theory. The holographic dictionary between the theories relates the 
Yang-Mills and string coupling constants g^^i = 2TTg^. The temperature of the gauge theory 
is equivalent to the Hawking temperature of the black hole horizon: 

T=^. (2.16) 

In the supergravity approximation (corresponding to the limit ^ co^ g'^^^N^ ^ oo in 
the field theory), full information about functions li^{u},q) and Il^{u},q) can be obtained by 
solving the linearized Maxwell equations for the bulk electric field components E'r, El [29] 

Et+^jE't + ^^^,Et = 0, (2.17) 

-^ + 7^-^^^--0. (2.1B) 

where ' indicates a derivative with respect to x = 1 — rg/r^. We have also introduced the 
dimensionless quantities 

°=2if- <' = 2lr- 

An analysis of eqs. (2.17), (2.18) including the perturbative solution for small t», q can be 
found in [29]. Here we shall focus on a particular case of vanishing spatial momentum that 
admits analytic solution for arbitrary frequency. 

For q = 0, the components Et = El = E obey the same equation 



E" +^—E' + E = {). (2.20) 



f (1-^)/ 

Writing 

E{x) = (2 - x)-"'/2 F{x) , (2.21) 

where F{x) is by construction regular at the horizon x = 0, we obtain the equation 

„ 2iw + 2{x - 1) - (1 + i)wx + -x)- m((l + 2i) -x)) 

{x-2)x 2x{x-l){x-2) • V • J 

Two linearly independent solutions of eq. (2.22) are written in terms of the Gauss hypergeo- 
metric function 

. , , , (i+i)vo ( (1 + i)t» (l+i)t» X \ , ^ 

Fi{x) = (1 - x) — 2F1 i^l - , ; 1 - i»; ^^^^T) J ' (^.23) 

.TO.-, ^ (l-i)ro il-i)vo 



F2(x) = x*"* (1 - x)^2Fi (^1 - , ^ ' -A + ito; ^^ZTTj) " (^.24) 

To compute the retarded correlators, we need a solution obeying the incoming wave boundary 
condition at x = [15]. The correct solution is thus given by eq. (2.23). 
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The retarded correlation functions can be computed from the boundary supergravity 
action using the Lorentzian AdS/CFT prescription [15]. For vanishing spatial momentum, 
the result reads [29, 13] 

N^T^ ,. E'{x) 



n(a;) = 



lim 



*i E{x) 

Substituting the solution (2.23) into eq. (2.25) we obtain 



n(u;) 



zt» + to 



1 - i)K> 



(1 + i)t» 



(2.25) 



(2.26) 



where 'tp{z) is the logarithmic derivative of the gamma-function. The spectral function is 
given by 



2rp2 



NfT 



Im < i t» + t» 



(1 



z to 



(1 + i)to 



(2.27) 



Using the property of the digamma function il){z) — i^{—z) = — TrcotTr^; — 1/ the spectral 
function (2.27) can be written in a more compact form 



^l(a;) 



N'^T'^ TTto^ sinh Trto 



cosh 7rt» — cos 7rt» 



(2.28) 



These analytic results for the retarded Green's function (2.26) and the spectral function (2.27), 
(2.28) are new and allow their various features to be easily examined. The asymptotics of 
the spectral function for large and small frequency can be easily computed 



= " ] (1 + 2e 



-Trto 



^(w) 



4 



cos 7rt» 



+•••), 



oo , 



1 + 



6 



+ 



to 



0. 



(2.29) 
(2.30) 



As expected, the high frequency asymptotic coincides with the zero-temperature result for 
the spectral function [15] 

The retarded correlator (2.26) is a meromorphic function of to with poles located at^ 

vo = ±n-in, n=l,2, .... (2.32) 

The position of the poles coincides with the quasinormal spectrum of the fluctuations E{x) 
determined by the Dirichlet condition E(l) = 0. For each pole, the imaginary part has the 
same magnitude as the real one, and thus none of the singularities can be given a 'quasipar- 

ticle' interpretation. Indeed, as shown in fig. 3a, the spectral function is quite featureless, 
although not monotonic: its finitc-tcmpcraturc part, 'K{uj) — "K^^^ioj)^ exhibits damped os- 
cillations reflecting the diminishing influence of the sequence of poles receding farther and 
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(a) (b) 

Figure 3: The A/" = 4 SYM i?-current spectral function at zero spatial momentum ■K(a') (a) and its 
finite temperature part — N'^u'^ /IQn (b) in units of N'^T'^ /2 as a function of i» = uj/2'kT. 

farther away from the real axis in the complex frequency plane (fig. 3b). The oscillatory 
behaviour of the finite-temperature part of the spectral function is evident from eq. (2.29). 

Using the Green-Kubo formula (2.11) and the low frequency limit (2.30) of the spectral 
function at zero spatial momentum one finds the product of the i?-charge diffusion constant 
and the charge susceptibility 

«==^. (2.33) 

The susceptibility is determined from thermodynamics according to eq. (2.12). The depen- 
dence of the charge density on the chemical potential was found in [34]. For small //, one 
has 

n(/x) = ^^/i + --- , (2.34) 

and thus from eq. (2.12), H = We conclude that the R charge diffusion constant 

is given hy D = 1/2ttT, in agreement with the result of an earlier calculation [16], where 
the value of D was determined from the hydrodynamic pole of the longitudinal part of the 
correlator at small but nonvanishing spatial momentum. 

3. Adding flavour: D7-brane embedding and thermodynamics 

All fields in the = 4 SYM theory transform in the adjoint representation of the SU{N^) 
gauge group. Fields transforming in the fundamental representation of the gauge group can 
be introduced in the gravity dual by inserting a second set of D-branes in the supergravity 
background [5]. In particular, we consider the decoupling limit of the intersection of black 
D3-branes and Nf D7-branes as described by the array: 

0123456789 
L»3 X X X X (3.1) 
D7 xxxxxxxx 



^The exact location of the poles was previously found in [33] using the continued fraction method. 
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The dual field theory is now an = 2 gauge theory consisting of the original SYM theory 
coupled to Nf fundamental hyper multiplets. Taking the decoupling limit with N( <^ N^, the 
D7-branes may be treated as probes in the black D3-brane geometry (2.15). The holographic 
framework has been used extensively to study the thermal properties of the M = 2 gauge 
theory at large A^^ [8, 9, 7, 10]. In particular, it was found that the fundamental matter 
undergoes a phase transition characterised by the dissociation of the mesonic bound states, 
as will be reviewed below. 

Following [7, 10], it is helpful to introduce a dimensionless radial coordinate p, related to 
the coordinate r via 

{top? -- 



+ \A 



In this case, the metric for the black D3-brane geometry (2.15) becomes 

L2 



f 



+ 



where f{p) = 1 — f{p) = 1 + 1/p'^ and 



dnl = de^ + sin^ ed^s + cos^ e d(j)'^ 



(3.2) 



(3.3) 



(3.4) 



The worldvolume directions of the D3-branes by the coordinates {i, x*}. The probe D7-branes 
fill these coordinates, as well as wrapping the in (3.4) and extending in the radial direction 
p. In general then, the D7-brane embedding is specified by giving its profile in the angular 
directions, 6 and (j). For simplicity, we fix = for our background embeddings. Then 
requiring translational symmetry in the 0123-space and rotational symmetry in the 4567- 
directions motivates us to consider 6 = 9{p). In fact, it is more convenient to work with 
x{p) = cos9{p) rather than 0. The induced metric on the D7-branes is then 



, 2 1 (ropy 

= 2 y-T) 



f2 

—dt'^ + fdx^ 
f 



I? 
p2 



2-2 



+ 1^+1 



dp'' + L\l-x^)d^l, 



(3.5) 



where % = dx/dp. The D7-brane action follows as 



where the normalisation constant \sM = NfTj^Yro^n^/AT with = 27r^ and 
The equation of motion for x{p) is then 



(3.6) 

D7 - 27r/(27r4)^5. 



do 



1 - 



1 



p5(l 



+ P' I 



P 



which implies that the field x asymptotically approaches zero as 

m c 
P P'^ 



l\ 3x(l - x') + 2p\x^ 

,8 



X" 







(3.7) 



(3.8) 
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The operator dual to x is the supersymmetric extension of the quark mass term, defined in 
(A. 2). Holography then relates the dimensionless constants m and c to the quark mass and 
condensate via (A.8) and (A.9). Eq. (A. 8) implies the relationship m = M jT between the 
dimensionless quantity m, the temperature T and the mass scale 

Here Mg^p is the meson mass gap in the D3/D7 brane theory at zero temperature [11]. 

In the limits of large and small m it is possible to find approximate analytic solutions 
for the embeddings - see [10]. However, for arbitrary m we numerically integrated (3.7) - 
see [7, 10]. In the present case, we are studying the high temperature phase and so we arc 
interesting in the black hole embeddings, which are found by imposing the following boundary 
conditions on the event horizon p^^^ = 1: x — Xo *iid dx/dp = ioi < Xo < ^- In order 
to compute the constants m, c corresponding to each value of xo, we fitted the numerical 
solutions to the asymptotic form (3.8). 

3.1 Thermodynamics of the brane 

As described above, introducing D7-brane probes into the black D3-brane background cor- 
responds to adding dynamical quarks to the gauge theory. The resulting theory has a rich 
spectrum of quark-antiquark bound states or mesons. As these mesons are dual to strings 
with both ends on D7-branes, the mesons can be studied by examining the fluctuations of 
the D7-branes, e.g., the spectrum of the lowest-lying mesons can be found by computing the 
spectrum of fluctuations of the worldvolume fields on the D7-branes - see e.g., [11, 12, 35]. 
At finite temperature, there are two classes of embeddings for the D7-branes and a first-order 
phase transition that goes between these classes. For temperatures below the phase transi- 
tion T < Tfund, the D7-branes close off above the black hole horizon (Minkowski embeddings), 
while above the transition (T > Tf^^^), the D7-branes extend through the horizon (black hole 
embeddings). In the gauge theory, the most striking feature of this transition is the change in 
the meson spectrum [7, 10]. Below the phase transition, the spectrum of mesons has a mass 
gap and is discrete while above, the mesonic excitations are unstable and are characterised 
by a discrete spectrum of quasinormal modes [36] . 

The above phase transition has been studied in detail in [8, 9, 7, 10] and we review a few 
salient facts here needed for later discussions. The Minkowski and black hole embeddings are 
separated by a critical solution that just touches the black hole horizon, as depicted in fig. 1. 
Plotting the quark condensate c as a function oim = M /T reveals that c is not a single- valued 
function of m and that the two families of solutions spiral around the critical solution - see 
fig. 4 of [10]. The physical solution corresponds to the embedding which minimizes the free 
energy density of the D7-branes. 

A plot of the free energy F as a function of temperature near the phase transition is 
given in fig. 4 and this shows the 'swallow tail' which is typical of first order phase transi- 
tions. Starting from low temperatures, we follow the blue dotted line depicting Minkowski 
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embeddings to the point where this Une intersects the soUd red hne for black hole embed- 
dings. The phase transition occurs at this point, and the physical embedding jumps from a 
Minkowski embedding, a finite distance from the black hole horizon, to the black hole embed- 
ding with xo ~ 0.94. At this temperature, the quark condensate, entropy, and energy density 
each exhibit a finite discontinuity, indicating that the phase transition is first order. 
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Figure 4: The free energy F/MT versus temperature T/M for a D7-brane in the black D3-brane 
geometry where N = XNfNJT^ /'i2. The blue dashed (red continuous) curves correspond to the 
Minkowski (black hole) embeddings. The values of xo for certain black hole D7-brane embeddings 
are noted for future reference. The phase transition is indicated by the vertical dashed line labelled 
Xo = 0.94. 

It is interesting to ask whether the D7-brane embeddings beyond the phase transition, 
e.g., between Ai and A2 on the black hole branch, still represent metastable configurations. If 
so, the corresponding states of the gauge theory might be accessed by a process analogous to 
'super-cooling' the system. Examining the specific heat Cv = dE/dT reveals that Cv becomes 
negative as the curves {e.g., of the condensate, entropy or energy density as a function of 
T) spiral around the critical solution, indicating that the system should be unstable for 
these embeddings. In particular, the specific heat first becomes negative at A2 on the black 
hole branch and ^3 on the Minkowski branch. Examining the scalar fluctuation spectrum 
of the D7-brane Minkowski embeddings (corresponding to the meson spectrum in the low 
temperature phase of the dual gauge theory) reveals that a dynamical instability appears 
precisely at ^3: At the first kink in the free energy, the lowest-lying scalar mode on the D7- 
branes becomes tachyonic. In fact, at the second kink, the second lowest-lying scalar mode 
becomes tachyonic and new tachyonic modes seem to appear at each such kink [10]. 

Hence the behaviour on the Minkowski branch is clear: Continuing along the Minkowski 
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branch past the phase transition, the system exhibits a dynamical instabihty, which matches 
the thermodynamic prediction, at the point in fig. 4, which is the first kink in the free 
energy on the Minkowski branch and corresponds to the first turn-around in the spiral. Hence 
while these configurations remain metastable between Ai and A^, all of the embeddings 
beyond ^3 are simply unstable. In fact, more and more instabilities appear as the embeddings 
approach the critical solution, as described above. 

On the black hole branch we expect similar phenomena. At point A2 in fig. 4, the specific 
heat becomes negative, indicating a thermodynamic instability. Though the full calculation of 
the quasinormal spectrum remains to be performed - preliminary results appear in [36] - we 
will see in section 4.2.2 that the scalar spectral function provides evidence that new tachyonic 
modes again appear at each turn in the spiral along the black hole branch. Appendix D 
presents a complementary analysis which also supports the appearance of tachyons in the 
quasinormal spectrum of the scalar fiuctuations. Hence we expect that only the configurations 
between Ai and A2 on the black hole branch represent metastable states of the gauge theory. 

4. Spectral functions for excitations of fundamental fields 

In this section, we compute spectral functions for excitations of fundamental fields in the high 
temperature phase of the theory, T > Tfu^d, by studying vector and scalar fluctuations of the 
D7-brane probes. The details of the holographic dictionary relating the fiuctuations of the 
probe branes to the hypermultiplet operators of the gauge theory arc described in appendix 
A. In each case, we begin by considering modes that are constant on the internal 5'^ but 
then extend the analysis to include modes with nontrivial angular momentum on this space. 
The latter modes are dual to higher dimension operators which are 'charged' under the global 
S0(4) symmetry, as outlined in appendix A. 

4.1 Vector 

In the gravity dual, the vector is one of several possible excitations of the worldvolume gauge 
field on the D7-brancs. These modes arc characterised as having only At^x,y,z nonzero with 
Ap = Ags = [12]. The holographic dictionary outlined in appendix A reveals that the vector 
is dual to the current which is the conserved current corresponding to the diagonal U{l)q 
of the global fiavour symmetry. 

The full action for the gauge fields on a D7-brane contains the Dirac-Born-Infeld (DBI) 
action plus a Wess-Zumino term, however, for gauge fields with Ap = A^s = only the DBI 
portion of the action is relevant. Further since we only study linearized fluctuations about 
the background, the gauge fleld action is only needed to quadratic order, which is simply 



where the Latin indices run over the D7 worldvolume directions and gab is the induced metric 
on the D7-brane given in (3.3). 




(4.1) 
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Assuming that the gauge field is independent of the coordinates on the , we can easily 
reduce (4.1) to an effective action in five-dimensions. The induced metric in these directions 
is 



ds\g) 



2 \ L J 



-ijdt^ + fdxl 



+ — 



and so the determinant of the full induced metric (3.3) can be written as 



9lff{p) 



L3(l-x^)f. 



(4.2) 



(4.3) 



Here is the determinant of the metric on the of unit radius and g^ff is a radially- 
dependent 'effective coupling'. Integrating over the three-sphere, the action (4.1) reduces 
to 



S = -- 



I 



Vl'iT^^Nt / dtd•'xdp^/~g 



a 13 
?e//(P) 



where a,(3 = t,x,y,z,p. Of course, Maxwell's equations follow as 



(4.4) 



da 



9eff 



0. 



(4.5) 



Using the equation of motion (4.5), the action (4.4) can be written as 



(27r< 



|2\2 



-TorQ^Nt / dx^dpd, 



(4.6) 



Retaining only the terms at the p-boundaries and using the metric (4.2), this becomes 



(2<^)^ ^ ^ ro\. -4 



fp%i-x') 



2\2 



, AidpAi - —AtdgAt 



, (4.7) 



where i is summed over x,y,z. Following [29], we take the Fourier transform of the gauge 
field and with = (— q, 0, 0), 



A, 



^e-^'+''^^A,{k,p) 



(4.8) 



(with Ap = 0, as discussed earlier), and the boundary action can be written as 

NfN^T^ f diodq 



24 



/dujc 



, A,{p,-k)dpAi{p,k) - J-At{p,-k)dpAt{p,k) 



We construct gauge-invariant components of the electric field: Ex = qAt + ujAx and Ey^^ = 
u)Ay^z- Note that in the language of section 2, E^ corresponds to the longitudinal electric 
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field El while Ey^z correspond to the transverse electric field Et- With these gauge- invariant 
fields, the action can be rewritten as (using eq. (4.15) below) 



NjN.T^ f diodq 



24 



(27r) 



/p'(l - X 



2\2 



Ex{p,-k)dpEx{p, k) 



-\ {Ey{p, -k)dpEyip, k) + E,{p, -k)dpE,{p, k)) 



(4.9) 



Focusing on the longitudinal electric field, we write 

Ekip) 



E^{k,p)=Eo{k) 



EkiPo 



(4.10) 



where it is understood that eventually the limit ~^ will be taken. We can then define 
the flux factor for E^ as [15]: 



NfNrT'^ 



24 



E_k{p)dpEk{p) 



(4.11) 



The usual AdS/CFT prescription tells us to evaluate it at the boundary p — oo to find the 
retarded Green's function for E^ [15]: 



G = -IT = 



NfNcT^ 
23 



2\2 



E^k{p)dpEk{p) 



Vl^F + TF - q^P/P)E.k{Poo)Ek{Po 



■ (4.12) 



The retarded Green's function for is the above expression times a;^, which for g = gives 



NfNcT^ 



Ekip) J 



(4.13) 



p— >oo 



upon using the asymptotic expansion (3.8) for x- Of course, this is the analogue of the 
expression in eq. (2.25) for the example discussed in section 2.3. The spectral function for 
g = is then 



^a:a;(w,0) = -2ImGxx(w,0) 



NiNcT"^ 



Im 



j^ dpEkjp) 
Ekip) 



(4.14) 



p— »oo 



In order to evaluate the spectral function, we must solve the equations of motion (4.5). 
For Ap = As3 =0 and A^ an s-wave on the S^, the equations for At and A^ are 



- g^'^'qA^ = , 



-9 ~tt, 



^eff 



-9- 



2 -r^r^'A, 



-9 ~tt~ 



{uqAt + u''A,)={}. 



(4.15) 
(4.16) 

(4.17) 
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Given the longitudinal field Ex = qAt + uiA^ as above, the system of equations (4.15)-(4.17) 
yields 



Ex + 



+ 5. In 



Substituting for the induced metric in (4.18), the equation of motion for E^ is: 
4wV7 



Ex = 0. (4.18) 



Ex + 



+8 



+ ^ 



//^^2 .2^ 



Ea: 



(4.19) 



Here we use the dimensionless frequency to and momentum q arc defined in cq. (2.19). 
Returning to the transverse electric field Et = Ey^^ j the equation of motion is 



-9 



2 -9^^<7™5pi?T 

9eff 



-9 



^ -5W((jV* + 9V'^)^r = 0. 
9h 



(4.20) 



For vanishing spatial momentum q = 0, this equation and that for Ex (eq. (4.19)) coincide. 
Thus, as discussed in section 2, for q = 0, the spectral functions are identical i.e., T^xxi^, 0) = 
'&yy{LO,0) = '&zz{^,0) and we denote these by 'H(a;) henceforth. 

We proceed to compute the spectral function ■^(a;) by solving the equation of motion 
(4.19) with q = 0. First we note that the case of massless quarks corresponds to the equatorial 
embedding of the D7-branes for which xip) = 0. Hence 5^^^ in eq. (4.4) is constant and the 
induced metric (4.2) matches precisely that of the AdSs black hole. Hence except for an 
overall normalization, the calculation of 11 (a;) is identical to that in the example discussed in 
section 2.3 and so in this case, it is possible to solve (4.19) exactly. We leave this exercise 
for the following subsection as it is a special case of the general analysis of charged vector 
operators, for which the case Mq = is also exactly soluble. 

For massive quarks (m 7^ 0), the embedding equation (3.7) must be solved numerically 
and hence it was necessary to numerically integrate (4.19) to solve for Ex. Near the horizon 
(p — ^ 1), eq. (4.19) implies that Ex ~ (p— 1)^*'^. Choosing the negative sign enforces incoming 
wave boundary conditions at the horizon. Thus, for each choice of quark mass, wc solved 
(4.19) numerically, taking Ex{p) = {p — l)^'^^F{p) and where F{p) is regular at the horizon 
with F(l) = 1 and dpF{l) = iw/2 for real w. 

The spectral function was evaluated using the numerical solutions for Ex and eq. (4.14). 
In the high frequency limit, the spectral function asymptotes to NfN^uj'^ /Air - see appendix 
C. Figure 5 provides plots of the finite-temperature part of the spectral function,^ 11 (a;) — 



*In using the wording 'finite temperature part of the spectral function,' we are adopting the language used 
previously for N = SYM in section 2.3 and ref. [29]. In the present case, this refers to the spectral function 
minus its high frequency asymptotics. 
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Figure 5: The finite temperature part of tlie vector spectral function, i.e., T^ — NtN^LiP'/ATi, in units of 
NfN^T"^ / 4:, versus vo — ijj/2-kT for various values of xo corresponding to different values of m = M /T . 
The upper plot shows values of xo corresponding to temperatures above the phase transition while the 
lower plot is for values of xo past the transition. The vertical dotted lines represent the mass of the 
lowest and first excited vector mesons in the low temperature (Minkowski) phase for a near-critical 
Minkowski embedding. 

NfN^u)'^ / Ait, for various D7-brane embeddings, specified by xo = xip = 1) (or equivalently by 
m). The upper plot shows the finite temperature part of the spectral function for temperatures 
above the phase transition: xo = ("t- = 0), xo = 0.1 (m = 0.1667), xo = 0.5 (m = 0.8080), 
Xo = 0.8 {m = 1.2026), xo = 0.94 (m = 1.3059) - the last of these corresponds to T/M^ 
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for the phase transition. Note that the x = Xo = 0-1 lines are virtuahy coincident. 
The lower plot shows the finite temperature part of the spectral function for values of xo 
corresponding to black hole embeddings after the phase transition, i.e., along the lines Ai to 
A2 and A2 to ^3 on the black hole branch in fig. 4. Note that as xo approaches 1, the finite 
temperature part of the spectral function displays high peaks. It is interesting to note that 
as Xo ~* 1) as well as growing sharper, the peaks in the spectral function become more closely 
spaced and move towards lower frequencies. For example, the peaks in the xo = 0.9999 line 
are much more closely spaced than those in the xo = 0-99 line. 

It is interesting to compare the positions of these peaks to the masses of the lowest vector 
mesons in the Minkowski phase [37]. The vertical dotted line at t» ~ 0.776 represents the mass 
of the lowest vector meson for a Minkowski embedding very close to the critical embedding. It 
seems that the position of the first peak of the spectral function is converging to a very similar 
value as Xo ^ 1 - certainly for xo = 0.9999, the first peak is very close to t» ~ 0.776. The 
second vertical dotted line at t» ~ 0.857 represents the mass of the first excited vector meson 
(n = 1,£ = 0) as the Minkowski embedding approaches the critical solution. In this case, it 
is likely that the second peak in the spectral function is approaching this value, but certainly 
it is not converging on this position as rapidly as the first peak. From the supergravity 
perspective, it is natural that these peaks are converging on the Minkowski phase spectrum, 
as described above, because near the critical embedding both the Minkowski and black hole 
embeddings of the D7-brane will be nearly identical except for very near the event horizon. 
The small horizon that still appears in the induced geometry of the near-critical black hole 
embeddings ensures, however, that there is a small imaginary frequency component of the 
quasinormal fluctuations. 

Note that in the Minkowski phase, t» ~ 0.98 corresponds to the mass of the lowest vector 
meson just at the phase transition. This is significantly above both masses quoted above 
near the critical solution - recall that in general the meson masses decreased as the critical 
embedding was approached [10]. Further then, this mass does not seem to be correlated 
with the positions of the spectral peaks for the black hole embedding at the phase transition, 
beyond being in the same general range. 

The high peaks in the spectral function for xo ^ 1 may be interpreted in terms of 
quasiparticle states because their width F and is much less than their frequency F ^ J7. 
Appendix D presents a complementary discussion which reaches the same conclusion. The 
appendix only explicitly discusses the pseudoscalar channel but the results for the vector 
are almost identical. In particular, the effective potential shown in appendix C develops a 
finite barrier at intermediate values of the radius as xo — > 1- This suggests the existence 
of metastable states in the corresponding Schroedinger problem which, as discussed in the 
appendix, would correspond to quasinormal frequencies with F <C in this regime. 

4.1.1 Charged vectors 

The J\f = 2 gauge theory under study here has an internal 50(4) = SU{2) x SU{2) global 
symmetry, which is dual to rotations on the D7-brane's internal S^. The vector modes 
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which are considered above are aU singlets under this symmetry. However, these operators 
only correspond to the lowest dimension operators in an infinite family of vector operators 
transforming in the {i/2,i/2) representation of the internal symmetry [11]. As outlined in 
appendix A, these operators are built up by combining the adjoint hypermultiplet fields 
(scalars) with the fundamental fields appearing in the singlet operators. 

Evaluating the spectral function for these vectors with £ 7^ follows closely the analysis 
in the previous section and so we only present the salient steps here. Of course, the first step 
is to consider an expansion of the world- volume vector in terms of spherical harmonics on the 
S\ 

A, = Y,y\S^)Al{p,x^), (4.21) 



with 



+ 2)y 



(4.22) 



where V^yj is the Laplacian on the unit three-sphere.^ Examining the eight-dimensional 
Maxwell equations arising from eq. (4.1), one finds that with £ 7^ (and T 7^ 0) cannot 
be set to zero in general. Hence the general analysis becomes somewhat more elaborate. 
However, if we focus on spatially independent fluctuations, i.e., the q ^ Q limit above, then 
both A^p ^ decouple and the calculations are greatly simplified. In this case then, the analog 
of eq. (4.9) becomes 



Si 



j 



2\2 



fp\i-x') 

^l-X^ + P^X^ 



El{p-k)dpEl{p,k) 



(4.23) 



where i is summed over x,y,z and Ef = coAj. Recall that with vanishing spatial momenta, 
there is no distinction between longitudinal and transverse electric fields in the language of 
section 2. 

Examining the asymptotic behaviour of any of the electric field components, we write 



E',{u;,p) = E',{u;) 



{ttTY , E,M 



2^2 ^~£;,,,(p. 



(4.24) 



where it is understood that eventually the limit poo ^ 00 will be taken. Note the factor of 
required to obtain the correct asymptotic behaviour - see appendix A. As above, taking 
variations of E^{k) then yields the the flux factor for Ef. This then leads to the following 
expression for the spectral function for Ai: 



^iicu) = -21mGiiLu) 



2e 



21+2 



2i+3 9pEe^M>) 



(4.25) 



p— »oo 



®Of course, the spherical harmonics for a given i are also labeled by two further SU{2) quantum num- 
bers, but we drop those as they are irrelevant in the following. Implicitly, our normalization is such that 

y'=\^=o{S^) = 1 and so Jd'nVh;y*J„,yl^ = 2n^6u' 5^m' Kn' ■ 
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with no sum on i. Instead for any value of i = x, y, z, the spectral functions are identical, i.e., 
'^xxi'^) — '^yyi'^) ~ '^izi^)^ because we are limiting our analysis here to the case of vanishing 
spatial momentum. Hence in the following, we denote these spectral functions by T^eiuj). 

In order to evaluate the spectral functions, we must solve the Maxwell equations arising 
from the eight-dimensional action (4.1). Expressed in terms of the electric field components, 
the relevant equation of motion is 



4/ 



+ 



dn 



X 



+ 



+ p'^x^ 



[i-x')f 



p2p 



- l{t + 2) 



(4.26) 



We now proceed to compute the spectral function (a;), first for massless quarks (m = 0) 
and then for quarks with a finite mass. 

Recall that the case of massless quarks corresponds to the equatorial embedding of the 
D7-branes for which x{p) = 0- In the previous section, we noted that for the £ = vector, 
the calculation of '&{u}) then becomes the same as that in our example of section 2.3. In 
particular, an analytic solution can be found for q = because it is possible to solve (4.19) 
exactly when x = 0. Here we show that in fact the general equation (4.26) for any i has an 
analytic solution in this case of massless quarks. 

Setting X = making the change of variables x = 1 — 2/ p^f = 1 — 2p^/(l + p^), the 
equation for the fluctuation Ei^^{x) is 

-2 i{l + 2) 



// /' / 

E'lo: + -J E'i^^ + 



Ei,oj = , 



(4.27) 



[{l-x)p 4(l-x)2/J 

where the prime denotes a derivative with respect to x. As in eq. (2.21), the solution is given 
bylO 

Ee,Ux) = x-^"'/^ (2 _ ^)-«>/2 ^(^) ^ (4 28) 

where the regular function F{x) is a straightforward generalization of the result (2.23): 



(l+i)tX) 

F{x) = {l-x)^^2Fi 



1 + 
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[l + i)to e {l + i)K> 



2 2 

The spectral function is then given by 



1 — ixx>: 



2t 



lim —^NtN.T 
e^O 2^ 



■2^+2 



Im/(x) 



1 + V7\ E{x,-m) 



1 



2(x - 1) 
E'{x,xo) 



Ed 



(4.29) 



(4.30) 



-t») E{1 

The right hand side of eq. (4.30) is independent of the radial coordinate [15] and thus can be 
computed at any value of x, e.g., at x = 0. We obtain 

2 



2% 



NfN,T'^^+^ sinhvrw 



1 + r 

2 2 2 



1 XX> IXX> 

1 _l \ 

2 2 2 



(4.31) 



^"Notc that near the horizon x ~ 2(p— 1)^ and so the small x behaviour here is consistent with the boundary 
condition at the horizon discussed for the numerical solution. 
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Eq. (4.31) shows that the poles of the retarded correlator corresponding to '^^(a;) are located 
at 

vo = ±(n + l + ^^ (iTi) , n = 0,l,.... (4.32) 

Note that there is an interesting degeneracy in the positions of these quasinormal modes in 
that their position only depends onn + £/2. This is reminiscent of the unexpected degeneracy 
found in [11], where the meson masses only depended on the combination n + £ {at T = 0). 
For £ = 0, eq. (4.31) reduces (up to the normalization) to the result (2.28). For odd and even 
£ > 0, respectively, eq. (4.31) can be written in the form 

^e=2n-i{u;) - -^N.N^T 27r[(2rz-l)!p cosh 7r«> + cos 7r«> ^ + (2^^^ 

2* [(2?^)!] cosh 7rt» — cos TTO) \ 4k* J 

where n = 1, 2, .... The asymptotics of the spectral function for large and small frequency are 

^^(o;) = lL^iV,iV^T2^+2 „,2£+2 + (_i)^2e-'^«' costtvo) (1 + O{l/vo^)) , (4.33) 

t» — oo , 

M^) = j^N,KT^'+^r\l + -)«,, «, ^ . (4.34) 
In particular, we have the £ = 1 spectral function: 

4 cosh TTW + cos 7rt» 

The large frequency asymptotics of T^i (u) is 

2 

^li(a;) ^ •Hi(a;) = yATfAT^T^ 2W (l - 26"'^"' costtto) + | , (4.36) 

where we have dropped 0{e^'^^) terms. Thus for sufficiently large values of to the finite 
temperature part of the spectral function, — ir^ NfN^T^vo'^ , exhibits damped oscillations 

around ir^NfN^T^/A - see fig. 6. Note that for £ > 2 the finite temperature part of the spectral 
function asymptotes t»^^+^ for large t» and thus the oscillatory behaviour again becomes a 
sub dominant effect. 

As in the previous section, for massive quarks (xo / 0), both the embedding equation 
(3.7) and the vector equation of motion (4.26) must be solved numerically. Solving for E£^^^ 
requires special attention to the boundary conditions near the horizon (p — > 1). As for 
the £ = case, the appropriate incoming wave conditions are imposed by taking Ei^^{p) = 
(p - l)-'^F{p) with F(l) = 1 and dpF{l) = ixo/2. 

The vector spectral function for £ = 1 is shown for various values of xo in fig- 7. For all 
values of xo; the £ = 1 spectral functions approach NfN^Lo"^ / IQir at large to - see appendix C. 



- 22 - 



1 2 3 4 5 

Figure 6: The finite temperature parts of tlie Afq = (xo — — 0,£— I vector spectral function 
(■Ki — A^fiVeW^'/lGTr) and of its high frequency asymptotics (4.36) ('Ki — 7VfA^cW'*/167r) (dashed blue line), 
in units of n'^NfN^T'^/S, versus t» = ui/2ttT. Note the figure also demonstrates the precise agreement 
between the numerical results (red dots) and the exact result (solid black line, which is essentially 
invisible above). 

While this common behaviour is not clear in fig. 7, it can be seen by going to larger i». Note 
that the spectral functions in the upper plot, which correspond to values of the M^/T above 
the phase transition, seem to be essentially featureless. In contrast, the lower plot shows that 
as the critical embedding is approached with Xo ~^ 1 some peaks are appearing in the spectral 
function. The masses of the lowest two ^ = 1 vector mesons in the low temperature phase 
for a near-critical Minkowski embedding have been included in this plot as well. While these 
lines lie close to the first peak for the xo = 0.9999 spectral function, the peaks do not seem 
to be converging to these positions nearly as rapidly as was seen for £ = 0. 

The peaks in the spectral function for Xo ^ 1 may again be interpreted in terms of 
quasiparticle states when their width T and is much less than their frequency 0: F ^ (7. 
Hence, as discussed for i = 0, it appears that the quasinormal frequencies are approaching the 
real axis in this regime. However, we stress that this approach is occurring much more slowly 
for the i > modes. In particular, the spectral function remains essentially structureless for 
Xo = 0.94, which corresponds to the phase transition between the black hole and Minkowski 
embeddings. Therefore the mesonic states corresponding to the higher-^ vector operators 
dissociate immediately at the phase transition. 

Note the complementary discussion in Appendix D would lead to similar conclusions. In 
particular, a barrier in the effective potential develops as in the 1 = analysis but one must go 
to values of xo much closer to one when ^ > 0. Hence metastable states in the corresponding 
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Figure 7: The vector spectral function for £ = 1 in units of tt^ N{N^T^ versus to. The upper plot 
shows values of xo corresponding to temperatures above the phase transition while the lower plot is 
for values of xo past the phase transition. In the lower plot we focus on values of v» for which the 
spectral function shows structure. As usual, the vertical dotted lines represent the mass of the lowest 
and first excited vector mesons for £ = 1 in the low temperature (Minkowski) phase for a near-critical 
Minkowski embedding. 

Schroedinger problem would only appear for xo ™- this regime very close to one. 
4.2 Scalars 

We now turn to scalar and pseudoscalar excitations of the fundamental fields. In the dual 
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gravity picture, these correspond to scalar fluctuations of the D7-brane probes in the black 
D3 geometry (3.3) about the fiducial embedding given by dv{p)' 



+ ScPia") 



(4.37) 



where a" denotes the D7-branes' worldvolume coordinates. Appendix A describes the holo- 
graphic dictionary relating the scalar 59 and the pseudoscalar 6(/) to the corresponding gauge 
theory operators. We present the analysis for general i modes but our results will focus on 
the spectral functions of lowest dimension operators with £ = 0. 

The pull-back of the bulk metric (3.3) to the D7 worldvolume is: 

2L'X 



where 



ds\g) - 



-.da{5e)dx''dp + l2 [da{5e)db{8e) + x^da{5(t))db{5(t))\ dx^dx^ (4.38) 



-Cdt^ + fdxl 



+ — 
P L 



1 + 



2 • 2 

P X 



dp^ + p^ sin^{ey + 5e)d^l 



(4.39) 

As before, we've put x{p) = cos 6y{p). Using the DBI action and retaining terms only to 
quadratic order in the fluctuations, the Lagrangian density is 



C = Co + - 
3 



X 



1-X')damdbi5e) 



(4.40) 



L 



2 1-;^2 + p2^2 ■ 2 



1 - + P^X^ 



+ x'da{Scl^)d,{5<l>) 



where is the metric (4.39) with 66 = 0, Cq is the Lagrangian density for the fiducial 
embedding % (given in equation (3.6)), and the boundary terms Ci and £,2 are 



Ci = dp 
= dp 



p'ffji-xy/'x ^ 

a/1 - x^ +P^X^ 
3 pV/(l - X')XX 
2 ^1 - x^ + P^X^ 



(4.41) 
(4.42) 



We eliminated terms linear in 66 by integrating by parts and using the equation of motion 
(3.7) for X- 

The equations of motion for the fluctuations follow from (4.40) as 



On 







for 6(j) and 



^3p-'ff{l-x')\ab 



a/1 -X^ + P^X^ 



9fdbiS6) 



+ 3^^^£mZ^66 = 
y/l - X'^ + P^X^ 



(4.43) 



(4.44) 



for 66. 
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4.2.1 Pseudoscalar 5^ 

The relevant portion of the action (4.40) for the pseudoscalar Scj) is 



J <fada [v^pV/(1 - xWl-X^^ + P^yf 5</'W] (4.45) 



where we've integrated by parts and used the equation of motion (4.43). To evaluate the 

spectral function we only need the complex part of (4.45) and hence in the following wc 
retain only the term involving the p derivative. Expanding the fluctuation in terms of spherical 
harmonics on the of unit radius, 



5^ = Y,y\s'')5Mp,x''), 



(4.46) 



the term needed to evaluate the spectral function for the ^th mode is 



I 



pV/(i-xW 



We take the Fourier transform of d(f)i with A; = {—u, q, 0, 0), 



(4.47) 



p— >oo 



diodq 
(2^' 



-iuit+iqx 



,k), 



and write 



6(pi{p,k) = d(p)l{k)^^p^ 



2^/2 '^^VeMPoo) 



(4.48) 



(4.49) 



where the limit poo — ^ oo will eventually be taken. Note the factor of p^ required to obtain 

the correct asymptotic behaviour 5(j)i{p^,k) = S(f)^{k)p^^ - see appendix A. We can then 
define the flux factor for the £th mode as 



2^+6 



Vl - + p'id 'Pe,-kiPoo)Ve,kiPoc) 



(4.50) 



/!>— >00 



The retarded Green's function is then G = —IT [15] from which we obtain the spectral 
function 'H = — 2Im G for (; = as 



^^,(u;,0) = -^AiVfiYeT^^+W Hm Im 



(4.51) 



where we have simplified using eq. (3.8). 

For the ^ = mode, the holographic dictionary in appendix A describes how the variation 
6(1)^ (k) introduced an insertion of the operator Mq O^p. Hence to get the spectral function for 
the dimension-three operator O^, we should normalize the spectral function by an extra factor 
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1/Mq. We expect that this should also hold for £ > 0, in which case the operator O^p^ has 
dimension £ + 3. Recalling that m? = 4M^/XT^, we arrive at: 




Using (4.48) and (4.49), the equation of motion (4.43) becomes 



pV/(i-xW 

a/1 - + p2^2 





£{£ + 2) Pf,fe = 0.(4.53) 



Near the horizon (p ^ 1) we impose incoming wave boundary conditions so that taking 



we were able to solve (4.53) numerically to evaluate the spectral function (4.52). The high 
frequency asymptotics of the spectral function are described in appendix C. 

Figure 8 provides plots of the finite temperature part of the spectral function, — 
NfN^uj'^ / Air, for the pseudoscalar dcp, £ = 0, for various values of xo- Qualitatively the results 
are the same as for the vector spectal function shown in figure 5.^^ The quasiparticle peaks 
in spectral function quickly dissipate above the phase transition, i.e., for xo < 0.94. High 
sharp peaks develop as xo — 1- As before, the position of these peaks may be compared with 
the masses of the lowest pseudoscalar mesons on the Minkowski branch. The vertical dotted 
lines mark the masses (to ~ 0.770 and 0.849) of the lowest two 5(j) mesons (with £ = 0) for a 
near-critical Minkowski embedding. Note that the first peak in the xo = 0.9999 line is nearly 
centred on the first value of t». The second peak of this spectral function also seems to be 
converging towards the mass of the next meson. Hence as in the vector channel, we see that 
the spectra of these pseudoscalar fluctuations in the Minkowski and black hole phases seem 
to converge as we approach the critical embedding. 

As before, the sharp peaks which develop in the spectral function as xo approaches 1 
may be interpreted in terms of quasiparticle states. Again the complementary discussion of 
the quasinormal spectrum, given in Appendix D, leads to the same conclusion. 

4.2.2 Scalar 69 

The derivation of the spectral function for the scalar 66 is entirely analogous to that for the 

^^In the interests of space wc do not include plots of the pseudoscalar spectral function for £ > here, 
however, the i = 1 plot closely ressembles that for the vector, shown in fig. 7. 



re,k{p) ^ (P - l)-^"* l + -{p-l) + 0{p - If for p ^ 1 , 



(4.54) 
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Figure 8: The finite temperature part of the Scj), i = 0, spectral function, "H^ — NfNcUp' /A-k, in units 
of NfN,T^/A versus t» = uj/2TrT for various values of xo- The upper plot shows the spectral function 
for values of xo corresponding to temperatures above the phase transition while the lower plot is for 
values of xo past the transition. The vertical dotted lines represent the masses of the lowest two 
pseudoscalar mesons for a near-critical Minkowski embedding. 



pseudoscalar. The portion of the action for the 



fluctuations is, from (4.40), 



(Fa 



^2 _^ p2-^2 



^2 _^ ^2^2 
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where we've integrated by parts and used the equation of motion (4.44). As discussed above, 
to evaluate the spectral function we only need the imaginary part of the Green's function and 
hence of this action. Thus, only the p derivative term from the second line is needed. We 
expand the scalar in terms of spherical harmonics on the (as in (4.46)), take the Fourier 
transform (as in (4.48)), and express the ^th mode as 



(4.55) 



where we will eventually take the limit poo —>■ oo. Note that the factor of p^^ is inserted to 
obtain the correct asymptotic behaviour - see appendix A. 

Following the same procedure as with the pseudoscalars, we identify 



2e 



2^+6 



XNtNX- 



p'ffil - X 



2\3 



p''Jr^ne^-k{p)dpne,k{p) 



(1_ ^2 + ^2^2)3/2 ne,^k{Poo)ne,k{Po 



(4.56) 



The spectral function then follows as 



^3+2e dp'^iMp) 
T^t,k{p) 



(4.57) 



where we've used (3.8) to simplify. 

Now recall % = cos^ and asymptotically x — f^/p where m is determined by gauge 
theory quantities in eq. (A.8). Note that asymptotically we can relate a variation in 9 with 
a variation in 5x = —59. Hence, a variation of the coefficient of the operator Om in the 
gauge theory action (i.e., figuratively we might say 5M^{k)) corresponds to {\/\T /2)69o{k) 
in cq. (4.55). In the correlator (4.57) two factors of 59q have been stripped off, so in order 
to normalize the correlator so that only the variations of the gauge theory coefficient are 
removed, we should multiply by a factor of 4/Ar^: 

r-2e. 



4 TT 

= ^ ^e(a^, 0) = 



,N,N^T'^^+'^ lim Im 



^3+2£ ^P^fc(p) 

Tlk{p) 



(4.58) 



With the Fourier transform of 59 and using the notation (4.55), the equation of motion 
(4.44) for 59 becomes 



do 



p'ff{l-X 



2\3 



(1_;^2 + ^2^2)3/2 -P'-^-^ 



+ 



p'//(i-x') 

a/1 -X^ + p2^2 



(1-x^) (f 



p'f \f 



7t» 



q2 -(£ + 3)(^-l) 



7^^,fe = 0.(4.59) 



As with the vector and pscudoscalar, we set q = and impose incoming wave boundary 
conditions at the horizon, requiring that the field behave as 



ne,k{p)-{p-l) 



-ItO 



ivo , 



l + :^(p-l) + 0(p-l)2 



(4.60) 
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Figure 9: The finite temperature part of tlie £ = scalar (69) spectral function, 'Kg — NfN^uj'^ / Att , in 
units of NiN^T^/A for xo < 0.94, corresponding to temperatures above the phase transition. 

near p = 1. 

We solved (4.59) numerically and evaluated the spectral function using (4.58). The high 
frequency asymptotics of the spectral function appear in appendix C. Plots of the finite 
temperature part of the s-wave spectral function, 'Sg — NfN^uj'^ / An are provided in fig. 9 
for D7-brane embeddings corresponding to temperatures above the phase transition. The 
spectral function shows no high peaks and little structure at temperatures above the phase 
transition. 

Figure 10 provides plots of the spectral function for values of 0.94 < xo < ^, corresponding 
to black hole embeddings past the phase transition, i.e., continuing along the black hole 
branch in fig. 4 past point Ai. For 0.94 < Xo < 0.96, prior to the first kink in the free energy 
(between Ai and A2 in fig. 4), no striking peaks appear in the spectral function. However, for 
Xo = 0.9621, point A2 in fig. 4, a very high peak appears in the spectral function, centred on 
uj = 0. Taking a value of xo slightly larger (smaller), say xo = 0.964 (xo = 0.96), fig. 10 shows 
that the peak is diminishing and is centred on a small but nonzero value of uj. A bit further 
away from the first kink, e.g., xo = 0.97,0.98 no peak is evident. Following the D7-brane 
embeddings to the second kink, which occurs for xo = 0.99973885, the same behaviour is 
evident: Near this value of xo a small peak starts to appear in the spectral function and at 
Xo = 0.99973885 a high peak, centred on uj = appears. As we will discuss in section 6, this 
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Figure 10: The finite temperature part of the scalar {69) spectral function for ^ = in units of 
for values of xo past the phase transition. The lower plot focusses on the region near t» = 
where a peak appears in the spectral for xo = 0.9621, corresponding to the first kink in the plot of 
the free energy versus temperature. 

behaviour is a result of quasinormal eigenfrequencies crossing the real axis from the lower to 
upper half of the complex w-plane. As a result, these black hole embeddings become unstable 
beyond xo = 0.9621, in precise agreement with the thermodynamic discussion of section 
3.1. Ref. [36] examines the quasinormal modes in this channel directly and find qualitative 
evidence of this behaviour, as well. 
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At first sight, tiie qualitative difference in the behaviour of the i = scalar spectral 
function from the previous cases may seem to be at odds with the expectation that the 
fluctuation spectra of Minkowski and black hole phases should converge as they approach the 
critical embedding. However, we must recall that it was precisely the i = scalar modes 
that also realized an infinite family of instabilities on the Minkowski branes [10]. Hence we 
expect that the spectra of the two phases are again converging but now on the imaginary 
frequency axis. Hence the absence of a series of quasiparticle peaks should in fact be the 
expected result. 

Note, however, that in the analysis of the Minkowski phase [10], the ^ > scalar modes 
remained stable. Hence we might expect to see the appearance of quasiparticles in the spectral 
functions for these channels on the black hole embeddings. Fig. 11 provides a plot of the scalar 
spectral function for the £ = 1 mode. The spectral function shows very little structure for 
any values of XO; or, equivalently, m. However, we certainly found no evidence for ^ = 1 
of quasinormal eigenfrequencies crossing the real axis from the lower to upper half of the 
complex w-plane. Instead from the general arguments above, we expect that if this £ = 1 
spectral function was studied more intensively that in fact quasiparticle peaks would appear 
very close to xo = 1- 

We close this section with one final observation. While the pseudoscalar equation of 
motion (4.53) is singular in the massless limit, i.e., xo 0; the spectral function should 
have a well defined limit. Further if comparing the X ^ ^ limits in the pseudoscalar and 
scalar channels in figures 8 and 9, we find that they converge on the same spectral function 
in this limit. The fact that these spectral functions coincide in this limit is a reflection of the 
restoration of an additional U{1) global symmetry, corresponding to rotations in the 89-plane 
in the array (3.1) when the quark mass vanishes. In the massless limit, this symmetry relates 
the two scalar operators. 

5. Diffusion constant for 'light' quairks 

The worldvolume gauge field is dual to a conserved current in the dual gauge theory. One then 
expects to see the diffusion of the conserved charge, i.e., quark charge, according to Fick's 
law with a certain diffusion constant. This expectation can be confirmed in a holographic 
context [17, 16] and in fact, the computation of the diffusion constant can be performed in a 
number of different ways. In the present D3/D7 brane system, wc have explicitly computed 
the diffusion constant in three different ways: (a) using the membrane paradigm [17]; (b) 
the Green-Kubo formula; and (c) the lowest quasinormal frequency in the diffusion channel. 
In this section, we describe these different computations and our results confirm the internal 
consistency of the holographic framework, in that we show these different methods all give 
the same result. 

5.1 Membrane paradigm method 

The computation of the diffusion constant via the membrane paradigm was discussed in [17] 
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0.5 1 1.5 2 

Figure 11: The scalar spectral function for ^ = 1 in units of tt'^ NfN^T* / 8 versus to. The upper plot 
is for values of xo corresponding to temperatures above the phase transition while the lower plot is 
for Xo psLSt the phase transition. Note in the lower plot that the lines for xo = 0.9621,0.99,0.9999 
roughly coincide with that for xo — 0.94 and that there is no structure suggesting the existence of 
quasiparticle states. 



where explicit formulae for various transport coefficients in terms of metric components for 
a wide class of metrics were derived. There, the authors considered perturbations of a black 
brane background and a formula for the diffusion constant (eq. (2.27) in [17], also quoted 
here in eq. (E.5)) resulted from a derivation of Pick's law. An analogous computation can be 
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performed for the D7-branes' vector field for black hole embeddings and it gives 



12 



D 



1 



Vhs gxx\/—gttgpp 

2(l-A-o)"= 



-gtt) gpp- 



dp 



/a/I - + P'^X^ 

/2p3(l_^2)2 



(5.1) 



where, in the first expression, the metric g is the induced metric on the D7-branes (3.5) and 
h is the determinant of the metric on the of unit radius. 



DT 



0. 15 
.125 
. 1 
0.075 

0. 05 
. 025 



4 

T/M 



Figure 12: The diffusion constant D times the temperature T versus temperature T/M for D7-brane 
probes in the blade D3-brane geometry. The dotted vertical line marks the temperature of the phase 
transition. 



Using the numerical solutions for the embedding we numerically integrated (5.1) to 
find DT. The results are plotted in figs. 12 and 13. Fig. 12 clearly shows that asymptotically 
at high temperatures, DT approaches l/27r. This coincides with the result for the diffusion 
constant of R-charges in = 4 SYM [16]. At a pragmatic level, this coincidence arises 
because both results are constructed from correlators of a Maxwell field in a AdSs black hole 
background. As the quark mass is increased, the induced geometry on the D7-brane deviates 
from that of the background geometry. Hence one finds a departure of DT away from 1 /27r 
as the ratio T/M decreases. Close to the phase transition, there is a rapid decrease and 
DT = 0.036 ~ 0.226/27r at the phase transition. If we continue following the black hole 
branch beyond the phase transition, DT continues to fall and it also becomes a multi-valued 

^■^Note that the same method can be apphed to compute the diffusion constant for the gauge theory corre- 
sponding to the supergravity configuration of a Dq-brane probe in the near-horizon black Dp-brane geometry 
- see appendix E. 



function of temperature, as shown in fig. 13. The latter simply reflects the fact that multiple 
embeddings can be found for a single temperature in the vicinity of the critical solution. 




T/M 



Figure 13: Plots of the diffusion constant D times the temperature T versus temperature T/M for 
D7-brane probes in the black D3-brane geometry, zooming in on the spiral behaviour for temperatures 
near the phase transition. 



5.2 Green-Kubo formula 

As discussed in section 2.1, the diffusion constant may also be computed using the Green-Kubo 
formula (2.11) which relates the product of the diffusion constant D and the susceptibility 
H to the slope of the vector spectral function = 0) for w — > 0: Dr. = \\m.^^Q'^{ijj) /2uj. 
The susceptibility is H = duq/ d^i\^=Q where Uq is the charge density and is the chemical 
potential, both for fundamental matter (quarks or their supersymmetric generalization). In 
order to compute the susceptibility, one must consider the D3/D7 brane system with a finite 
chemical potential [38]. The susceptibility can be computed directly from eq. (2.22) of [38] 



- '■■r- , (5.2) 

'/(I - x')VfHi - x'f + 8*] 

which applies for any black hole embeddings. From the appendix of that paper, jl and d are 
related to ^ and via 

Combining these definitions, we have 

dd 4 dn„ 

(5.4) 



which, interestingly, is independent of the 't Hooft coupling A. 

Note from eq. (5.2) that /2 = is equivalent to d = which means that we can calcu- 
late dfi/dd\^^Q from this equation and take the inverse for the desired derivative. Hence a 



straightforward calculation yields 



dd 



and evaluating at d = gives 
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Combining eqs. (5.4) and (5.6), our final result is 



" dp, 
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(5.7) 



Note that in limit of massless quarks (ie., % = 0), these expressions have a simple form 
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d=0 
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(5.8) 



Numerically evaluating the low frequency limit of the spectral function 1^ and the sus- 
ceptibility (using (5.7) and (5.6)), we computed the diffusion constant using (2.11), and the 
results confirm those displayed in figure 12. 

Further, if we combine (5.7) and (5.1), we are lead to write 



(1 - y2)3/2 

47r 



(5.9) 



which, in view of (2.11), provides a simple analytic expression for the low frequency (a; — ^ 0) 
limit of the vector spectral function: 



(1 - y2)3/2 



27r 



(5.10) 



In a related analysis, the electric conductivity a of this system with fundamental quarks was 
recently computed from Ohm's law using AdS/CFT techniques [39]. In the limit of vanishing 
quark density and small external electric field, the result of [39] reduces to our formula (5.9) 
above, as it should according to the generalized Einstein relation a/e^ = D'E. (Note that 
is set to one in the conventions of [39].) 

5.3 Lowest quasinormal frequency (in the diffusion channel) 

The final computation of the diffusion constant comes from examining the hydrodynamic 
dispersion relation, corresponding to the lowest quasinormal frequency. At small three- 
momentum q, the diffusion constant can be extracted from: u) = —iDq^ + 0{q'^) [17]. In 
principle, the calculation of the quasinormal mode spectrum from eq. (4.19) proceeds as 
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follows: For p oo, eq. (4.19) implies that Ex — A + Bp~^ for some constants A,B. Nor- 
malizable modes will be those with ^ = 0. Hence one method to determine the quasinormal 
frequencies is to use a two-dimensional shooting method, i.e., solving (4.19) numerically with 
incoming wave boundary conditions at the horizon and then tuning the (complex) frequency 
to find a solution behaving as asymptotically. For small q wc solved (4.19) for various 
m to determined the lowest quasinormal frequency and our results for the diffusion constant 
are identical to those plotted in fig. 12. 

6. Discussion 

In this paper, we used holography to investigate various aspects of the high temperature phase 
of an = 2 super- Yang-Mills theory with fundamental matter. The holographic description 
consists of probe D7-branes in the near-horizon background of D3-branes (in the limit of large- 
and large- A with fixed iVf). In the high temperature phase, the D7-branes extend through 
the event horizon of the AdSs black hole, which describes the theory at finite temperature. 
In [7, 10], this phase was denoted as the black hole branch since the metric induced on the 
worldvolumc of the D7-branes is itself a black hole. Even though the latter geometry no 
longer obeys Einstein's equations, the analysis of the hydrodynamic physics found previously 
for bulk fields, e.g., [40], is readily transferred to the worldvolume fields on the D7-brane. 
Hence we were able to examine the spectral function for various mesonic operators in section 

4, following [30, 18], and we calculated the diffusion constant for the quark charge in section 

5, adapting techniques from [17, 16]. 

As reviewed in section 3, the induced geometry (3.5) on the D7-brane is determined by 
first solving for the embedding profile from eq. (3.7) with appropriate asymptotic boundary 
conditions (3.8). Given the complexity of eq. (3.7), these geometries are only known in general 
from numerical integration. However, there is one particularly simple case, namely that of 
zero quark mass. In this case, the embedding is trivial, i.e., X = everywhere, and the 
induced geometry (3.5) is precisely that of (the direct product of) an AdSs black hole (with a 
constant 5^). Hence our analysis of sections 4.1 and 5 reduces to studying a Maxwell field in 
an AdSs black hole geometry and the results precisely match those found previously for bulk 
gauge fields. For example, the quark diffusion constant (with Mq = 0) matches precisely with 
the i?-charge diffusion constant calculated in [16, 17]. Further the new analytic expression of 
the vector spectral function presented in section 2.3 was extended to an analytic result for all 
of the higher-^ modes of the worldvolume vector. 

As the quark mass is increased away from zero (or T/M^ decreases to finite values), the 
induced black hole geometry on the D7-brane begins to deviate from that of the background. 
In particular, the main differences arise near the event horizon where x is largest. For example, 
eq. (3.5) shows that the size of the and hence the induced horizon area shrinks as xo grows. 
Hence the physical properties of the fundamental fields were seen to depart (dramatically in 
some cases) from the standard results with the growth of the quark mass. Recall that some 
of the most interesting behaviour appeared as xo ^ 1; ^-c, approaching the critical solution 
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for which the effective horizon area vanishes. Hence this behaviour can be seen as a precursor 
of the phase transition to the low temperature phase or the Minkowski branch, in which the 
D7-brane smoothly closes off above the event horizon. 

In the low temperature phase (and in the limit of large N^), the spectrum of mesons is 
characterized by a discrete set of stable states [11, 10] and the spectral function is a series of 
(^-function peaks, as illustrated in fig.2a. A derivation of the spectral function for the scalar 
meson at T = appears in appendix B. These mesonic states correspond to open string 
excitations which are essentially living at the minimum radius of the D7-brane. Since in the 
high temperature phase the D7-branes extends through the event horizon, these states are 
destabilized. In this phase, the spectrum can be characterized by a discrete set of quasinormal 
modes in the effective black hole metric induced on the D7-brane. The spectral functions 
calculated in section 4 reveal interesting information about this quasinormal spectrum. We 
focussed on three particular operators, which are bilinears of the fundamental fields - the 
details appear in appendix A - corresponding vector, pseudoscalar and scalar channels. 

The behaviours of the vector and pseudoscalar spectral functions are very similar, as 
can be seen in figures 5 and 8. The following physical picture emerges from these plots: At 
very high temperatures, the spectral function closely resembles that for a vector in TV = 4 
SYM. (Of course, as discussed above, the bulk and worldvolume vector results are identical 
for r/Mq oo.) In this regime, the spectral functions show essentially no structure and 
the eigenfrequencies of quasinormal modes must all be deep in the (lower) complex plane. 
As the temperature decreases (with fixed quark mass), both the real and imaginary parts of 
a given quasinormal frequency decrease but the formation of peaks in the spectral function 
suggests that the imaginary part decreases more rapidly. At temperatures just above the 
phase transition, there are vector and pseudoscalar quasiparticles. Continuing along the 
black hole branch to even lower temperatures beyond the phase transition {i.e., following 
the black hole line through Ai in fig. 4), peaks grow very sharp and even more prominent 
indicating that the quasinormal modes have Re(a;) ^ Im(a;) in this regime. 

Appendix D presents a complementary discussion which reaches the same conclusion. 
Plots of the effective potential for the pseudoscalar (and vector) excitations in appendix D 
show a finite potential barrier developing at intermediate values of the radius as xo — ^ 1- This 
suggests the existence of metastable states in the corresponding Schroedinger problem which, 
as discussed in the appendix, would correspond to a quasinormal frequency with F CI i.e., 
the eigenfrequency approaches the real axis in this regime. Of course, while this intuitive 
picture developed from the effective potential matches the behaviour of the spectral functions 
discussed above, it only gives a very schematic picture of the quasinormal spectrum. Hence it 
would be interesting to develop more detailed picture with a full calculation of the quasinormal 
modes [36]. 

If we examine the positions of the peaks in the spectral functions more closely as the black 
hole embedding approaches the critical solution, it appears that the real parts of the quasi- 
normal frequencies roughly match with the spectrum of the lowest (vector and pseudoscalar) 
mesons on a near-critical Minkowski embedding. Hence one notable feature of the spectral 
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functions is that as xo ^ 1) tlie peaks are becoming sharper but also more closely spaced 
and moving towards lower frequencies. For example, in figure 8, the peaks in the xo = 0.9999 
line are much more closely spaced than those in the xo = 0.99 line. This behaviour is similar 
to what is seen for the dcj) spectrum for near-critical Minkowski embcddings: For these near- 
critical cmbcddings, the tower of masses appears to be collapsing to the mass of the lowest 
meson - see fig. 7 in [10]. As the phase transition occurs well away from the critical solution 
{i.e., Xo = 0.94 versus 1), the positions of the spectral peaks are not closely matched with 
the corresponding meson spectrum for the Minkowski embedding at the phase transition. Of 
course, both spectra still characterized by the same general mass scale, as given in eq. (3.9). 

We also examined the spectral functions for £ > in vector and pseudoscalar channels. 
These modes of the worldvolume fields correspond to higher dimension operators in the field 
theory, which are charged under the internal symmetry group SO {4) = SU{2) x SU{2) - 
see the discussion in appendix A. For these modes, the results are qualitatively similar to 
those for ^ = as one approaches the critical embedding. However, the rate at which the 
quasinormal frequencies approach the real axis is much slower - that is, the spectral functions 
only develop pronounced peaks very close to xo = 1- In fa-ct, these peaks are already washed 
out at the phase transition, i.e., xo = 0.94. Hence the corresponding mesons with ^ > do 
not survive as quasiparticles through the phase transition. The analogous observation applies 
to the excited mesons with £ = and n > 1. Examining the spectral functions in figs. 5 and 
8, one finds that only the first peak remains pronounced at xo = 0.94. Hence it seems that 
only the ground state mesons (with n = = i) can be said to survive the phase transition as 
quasiparticles. However, even these resonances have disappeared in the quark-gluon plasma 
by Xo — 0.8 or T ~ l.lTf^^^, where Tfu^d is the temperature of the phase transition. 

As shown in figures 9 and 10, the behaviour of the scalar spectral function is qualitatively 
different from that found in the vector and pseudoscalar channels, shown in figs. 5 and 8. 
As before, at high temperatures, the spectral function shows no distinguished structure, 
indicating that the quasinormal eigenfrequencies are all deep in the (lower) complex plane. 
As the temperature decreases, a small peak develops near the origin although it is still not 
very prominent at the phase transition. However, continuing to the D7-branc cmbcddings 
on the black hole branch for temperatures below the phase transition, this single peak grows 
and becomes extremely sharp and centred at a; = 0, precisely at the first kink in the free 
energy, i.e., xo = 0.9621. Beyond this point, the peak decays and moves away from a; = 0. 
Our interpretation of this behaviour is that the lowest (pair) of quasinormal frequencies 
approaches the origin and actually crosses the real axis at the point A2 in fig. 4. Continuing 
beyond this point, this eigenfrequency moves into the upper half plane, where it actually 
corresponds to an unstable mode. 

This interpretation of the behaviour of the scalar spectral function is confirmed by the 
qualitative analysis of the corresponding quasinormal modes in appendix D. Examining the 
effective potential for the scalar excitations shows that a negative potential well develops 
and grows as xo ^ 1- As discussed in the appendix, when this well is large enough, it can 
support long-lived 'bound' states for which (the real part of) the effective energy is negative. 
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These modes arc distinguished since > and further F < 0, so that these bound states 
correspond to instabiUties of the D7-brane. The spike (or pole) at a; = in the scalar 
spectral function discussed above results from the formation of the first bound state where 
the eigenfrequency crosses the real axis. Rcf. [36] studied the quasinormal modes in this 
channel and developed a qualitative picture which is in agreement with our results. 

Recall that the thermodynamic discussion of section 3.1 predicted that the system should 
become unstable at the point A2 in fig. 4 because the specific heat of the black hole branch 
becomes negative there. Hence our analysis above is in precise agreement with this result and 
it shows that the instability corresponds to unstable 'quasinormal' modes appearing on the 
D7-branes. 

In fact, we found the scalar spectral function also displays a spike at = at the second 
kink in the free energy, i.e., xo = 0.99973885. Hence it appears that the second lowest 
quasinormal mode becomes unstable at this point. It is natural to conjecture then that each 
time the free energy turns, a new 'tachyonic' mode appears in the scalar spectrum. In fact, 
the spectrum of scalar mesons on the Minkowski branch was found to display precisely this 
behaviour [10]. 

Again the spectral functions of mesonic operators which we calculated exhibited a number 
of interesting features, which had a clear interpretation in terms of the spectrum of quasi- 
normal modes. It would, of course, be interesting to confirm these behaviours by a detailed 
investigation of the quasinormal modes, as was begun in [36]. In the present paper, the 
spectral functions were only calculated for zero spatial momentum for computational simplic- 
ity. So another natural extension of this work is to consider the behaviour at nonvanishing 
spatial momentum. In particular, the spectral functions for general time-like and light-like 
four-momentum can be used to calculate photon and dilepton production rates, respectively 
[13]. An analysis of these results for the present M = 2 gauge theory has been made recently 
[41]. 

The other main result of this paper was the calculation of the diffusion constant for 
the quark charge. We used a number of techniques developed for bulk black holes in the 
calculation of the i?-chargc diffusion constant in the = 4 theory: the membrane paradigm 
method [17], the Green-Kubo formula [16] (which for us relied on previous studies of this 
system at finite quark density [38]), and the lowest quasinormal frequency [17]. It is gratifying 
that as demanded by the internal consistency of the holographic framework, all three of these 
independent calculations yield the same results [42], which are shown in figure 12. 

At very high temperatures [i.e., T jM ^ 1), the diffusion constant approaches ItxD T = 1, 
which as discussed above matches the i?-charge diffusion constant for the = 4 super- Yang- 
Mills theory [17]. As TjM decreases, the product DT decreases. Intuitively, we might 
understand this result as the rate of diffusion decreasing at a fixed temperature when the 
quark mass is increased. In fact, the decrease is remarkably small at first, e.g., 2'kDT = 
0.9 at T/M ~ 1.5. However, figure 12 then shows a dramatic decline as we approach the 
phase transition at T/M = 0.7658. At the phase transition, 2ttD T ~ 0.226, but if we 
continue following the black hole branch it appears that D T continues to decrease and would 
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vanish at the critical embedding. The membrane paradigm approach to calculating Z), as 
described in section 5.1, provides a straightforward understanding of this vanishing from the 
bulk perspective. Examining the expression in eq. (5.1), we see that in the critical limit, the 
integral remains finite but the prefactor vanishes because \/—g 0. Hence the dominant 
effect which causes D to decrease is the reduction of the effective horizon area of the brane 
geometry (3.5) as we approach the critical embedding. That is, (1 — Xo) — > in advancing 
toward the critical embedding. Further, this area vanishes precisely at the critical embedding 
confirming that D should vanish there. 

We have considered the M = 2 gauge theory to have N{ flavours of quarks. Hence it 
is worth noting that the results for the quark diffusion constant in the present holographic 
framework are independent of both Nf and N^. Of course, the same independence of was 
seen with the i?-charge [17]. This must certainly arise because we are working in the limit of 
large A^^ and large A. 

There has also been a great deal of interest in the diffusion of 'heavy quarks' in holo- 
graphic theories recently [46, 43, 44] - see also [45] and the references therein. In the present 

context, this simply refers to the quark diffusion constant in the low temperature phase where 
the D7-branes only approach to some finite distance away from the black hole horizon. In 
the low temperature or Minkowski embedding phase, a quark is represented by a funda- 
mental string stretching between the D7-brane and the horizon. As such, a heavy quark 
is holographically represented by a macroscopic object (on a similar footing with the probe 
D7-branes) and classically this object will remain at rest {i.e., it does not 'diffuse'). It is only 
when semiclassical effects are taken into account that the heavy quarks diffuse through the 
appearance of Hawking radiation from the effective black hole metric induced on the string 
worldsheet [43]. This process should be contrasted with the diffusion process in the high 
temperature phase which we have been considering here. As stressed above, in this phase 
the induced metric on the D7-brane is a black hole metric and so if quark number is injected 
into the system, the holographic description of diffusion is simply the classical process of the 
corresponding worldvolume excitations falling towards the event horizon. Given these two 
disparate descriptions, it is not surprising that the diffusion constant takes a qualitatively and 
quantitatively different form in the two phases. In the low temperature phase, the diffusion 
constant is governed by the heavy quark result [46, 43]^^ 

2t: DT = ^/^. (6.1) 

On the other hand, after the transition to the high temperature phase, we found above that 
2tt DT = 0(1). Note that the holographic analysis applies for strong 't Hooft coupling (i.e., 
A ^> 1) and so these results show the expected suppression of quark diffusion in the low 
temperature phase. 

^^Note we present this result with the same normahsation for the 't Hooft couphng used throughout this 
paper. Ref. [46] uses a convention such that A = 2A. This difference arises from the impUcit normahsation of 
the U{Nc) generators: Tr{TaTb) = dSat- The standard field theory convention used in [46] is d = 1/2 while 
our choice is d = 1, as is prevalent in the D-brane literature. 
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We comment on a possible puzzle with the above description of the diffusion of heavy 
quarks as due to semiclassical Hawking radiation. As such, the diffusion constant would be 
expected to vanish in the limit h ^ for the bulk theory. Now the standard AdS/CFT 
dictionary would associate Planck's constant h in the bulk supergravity with in the dual 
gauge theory [4]. However, above, we sec the result (6.1) is seen to be independent of A^e and 
D certainly does not vanish in the limit A,. — > oo. The resolution of this puzzle is that we 
have misidentified the correct 'semiclassical' nature of the diffusion process here. Above we 
observed that the heavy quarks diffuse because of the appearance of Hawking radiation in the 
effective field theory on the string worldsheet dual to such heavy quark. That is, the 'fields' 
on the worldsheet are the transverse coordinates describing the embedding of the string and 
so fluctuations in these fields (arising from Hawking radiation) corresponds to fluctuations in 
the position of the quark {i.e., diffusion of the quark). Now as usual, h for the worldsheet 
theory is identified with the inverse string tension, a' = More correctly, the dimensionless 
h in the nonlinear sigma model one the worldsheet is identified with the ratio of a' and the 
background curvature scale, i.e., hyjs — ^t/L'^ — I/VA- Now we see this physical picture 
matches precisely with the calculated result (6.1) and the limit A — oo is the semiclassical 
limit on the string worldsheet. 

The results for the diffusion constant in both the high and low temperature phases are 
combined in figure 14, where 2ttDT is shown as a function of A (for fixed MJT). The thick 
black curve shows a canonical result for the M = 2 gauge theory, which we are displaying 
for MJT = 200. The diffusion constant starts at very large values on the left in nearly 
perturbative results. The dashed 'perturbative' line is the extrapolation of the perturbative 
calculation of [47] for = 4 super- Yang-Mills at large N^. As the curve shows, we expect 
DT to make a transition to the A~^/^ behaviour of heavy quarks in the low temperature 
phase. If the quark mass was infinite, this behaviour would extend out to infinite 't Hooft 
coupling. However, for a finite mass, increasing A eventually takes the system to the high 
temperature phase. It would be interesting to understand the corrections to the heavy quark 
result (6.1), as we approach the phase transition. For M^/T = 200, the latter occurs at 
A ~ 9.4 X lO'' [7, 10]. At this first order phase transition, 2tt DT jumps discontimiously up 
to the 'light quark' curve at 27r DT = 0.226 and it quickly reaches the asymptotic value of 1 
as A continues to increase. 

Much of the recent interest in holographic calculations of the diffusion constant was 
generated by the possibility to compare these strong coupling results with experimental results 
for the QCD quark-gluon plasma measured at RHIC [48]. It is interesting then to place the 
RHIC results in the context of the phenomenology of the M = 2 gauge theory studied here. 
The bar labeled RHIC corresponds to 

Oistrong — 0-5 (or A — Stt) and the range 2tt DT — 3 
12, which is used in fitting the observed nuclear modification factor and elliptic flow amplitude 
for heavy (charm) quarks - see [48] for details. This value of the 't Hooft coupling is well 
away from the light quark behaviour of our canonical theory but also lies in an intermediate 
regime between the heavy quark and perturbative regimes. Hence a direct comparison of 
either approach for the N = 2 theory to the experimental data for QCD is difficult [47, 49]. 
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Figure 14: Sketch of 27r D T versus A for a canonical Af = 2 gauge tlieory witli AI^/T — 200. To tlie 
far left, we have the perturbative regime [47], which crosses over to the 'heavy' quark behaviour (6.1). 
Above the phase transition to the black hole phase at A ~ 9.4 x 10^, the curve reflects the 'light' quark 
behaviour of fig. 12. The light quark curves are also shown for M^/T = 2, 6, and 20. The horizontal 
dotted line marks the value oi2-K DT for the 'light' quarks at the phase transition. 



We emphasize that in this intermediate region om' canonical curve is simply an 'artistic' 
impression of the cross-over between these two regimes. 

For the canonical theory, we chose MJT = 200 so that the phase transition between 
the low and high temperature regimes took place for a value of A that we could confidently 
characterize as strong coupling. So we extend our discussion of RHIC results by noting that 
the charm quark has a mass of roughly 1500 MeV while the temperatures achieved in the 
RHIC collisions are in the regime 250 MeV. Hence in these experiments, we are considering 
Mq/r ~ 6. As illustrated in figure 14, the effect of reducing this ratio is to slide the light 
quark curve to the left. That is, the phase transition occurs at smaller and smaller values of 
A, e.g., Ac — 85 for MJT ~ 6. However, it seems that this critical value is still well away 
from the value of the coupling relevant for RHIC. Further it seems that the same will still be 
true for charm quarks even at the higher temperatures that might be achieved in the future 
at the LHC. Thus it is unlikely that a dramatic jump in the diffusion constant, such as that 
seen for our canonical theory, will appear in these experiments. It is noteworthy that in any 
event, the experimentally favoured values of the heavy quark diffusion constant are in fact 
above those calculated in the high temperature phase. Hence it would seem that a QCD 
phase transition would be characterized by a sudden decrease rather than a sudden increase 
in the diffusion constant. 

While our results for the diffusion constant do not seem relevant for 'heavy quarks', they 



-43- 



might be considered as that for 'hght quarks' in QCD. Note that our holographic model 
gave 2-ir DT < 1 where the effect of a finite quark mass was to give a slight (less than order 
one) reduction below the asymptotic value of 1. Appendix E extends the computation of the 
diffusion constant described in section 5.1 to more general holographic theories. In particular, 
figure 18 shows the results for the D4/D6-brane system, which is the basis for the construction 
of one holographic model which mimics QCD at large [50]. The results are similar to those 
above with 27r D T < 3/2 with finite mass effects giving a small reduction from the asymptotic 
value. Another interesting holographic model of a QCD-like theory comes from introducing 
D8 and anti-D8 probe branes in a D4-brane background [51]. The resulting diffusion constant 
is simply 2Tr DT = 1. In this model, the current quark mass is fixed vanish and so no finite 
mass effects appear. Wc might also recall the results for the i?-charge diffusion constants 
27r DT = 1 and 3/4, for the near-extremal D3- and D4-brane backgrounds [17]. Hence in 
all these cases then, we find that the calculations yield 27r DT = 0(1). This might suggest 
that in QCD, the diffusion constant associated with the light quarks falls dramatically at 
strong coupling, as compared to the perturbative results [52], but that they should saturate 
around this level. One might note then that these diffusion constants would be smaller than 
for heavy quarks at presently accessible energies but not much smaller. However, we must 
recall that these calculations are all performed at large A and large A^^ (with N(/N^ <C 1) and 
so it would of course be interesting to understand the corrections to these results at finite A 
and finite N^. 

In the present study, we focussed on calculating the diffusion constant for the overall quark 
charge. Of course, with Nf flavours of quarks with identical masses, the Af = 2 gauge theory 
has a global U{Nf) flavour symmetry and in the dual gravity description, the worldvolume 
theory of the D7-brane contains a nonabelian U (Af) gauge field. Our calculations have only 
considered the diagonal U{1) component of this gauge field. However, as noted in appendix 
A, one can easily examine the nonabelian SU{Nf) flavour currents with the corresponding 
components of the worldvolume gauge field — see, e.g., eq. (A. 6). Now in principle, we would 
describe the diffusion of the full set of flavour currents with a diffusion matrix Dab, rather 
than a single constant. However, our calculations of the diffusion constant in section 5 only 
relied on knowing the quadratic action for the dual gauge field and hence the nonabelian 
character of the gauge fields would play no role. Hence the diffusion matrix is, in fact, 
diagonal: D^b = DSab where D is precisely the constant determined for the U{1) charge. 
On general grounds, this is of course the expected result for the J\f = 2 gauge theory in the 
absence of any chemical potentials. With a nonvanishing chemical potential, the diffusion 
matrix will not take this simple form [53]. Similar comments to those above also apply for 
extending our calculations of the spectral function to operators that are no longer SU{N{) 
singlets. 
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A. Holographic dictionciry 

In sections 4.1 and 4.2, we focussed on three worldvolume fields: the gauge field and the 
embedding coordinates % and ^. We made use of the fact that the asymptotic behaviour 
of these D7-brane fields has a direct translation in terms of hypermultiplet operators in the 
gauge theory. In this appendix, wc elucidate this holographic dictionary in more detail. Our 
discussion provides more detail on the pscudoscalar operator since the other two operators 
have already been considered in some detail in [10, 38]. 

Let us remind the reader that a hypermultiplet consists of two Weyl fermions tp, ip and 
two complex scalars q, q. Of these, ip and q transform in the fundamental of the SU{N^) gauge 
group, while ■0 and q transform in the antifundamental. Further, with A^f flavours (of equal 
mass), the hypermultiplets transform under a global U {N{) ~ SU{Nf) x U (l)q symmetry. The 
charges of the fields under the diagonal t/(l)q are +1 for and q and -1 for and q. Hence 
the U{1)^ charge naturally counts the net number of quarks in a given state. Here and below, 
we follow the notation of [54]. 

The operators dual to A^, x and cf) can be determined by considering the interactions of 
the open strings on the D3/D7 array (3.1) before the decoupling limit [55], in analogy with 
the closed strings. Such an exercise leads to the following operators: 

Af, ^ = ^ai'i) + V^ct'^V''^ - i (q^Vq - {Vi'q)^) - i (q CD^q)^ - Vqq^^ , (A.l) 

X ^ 0^ = i'ipij + q{M^ + V2^^)q^ + q''{M^ + V2^^)q + h.c., (A.2) 
(f) ^ 0^ = ilj'ip + iV2q^^q^ + iV2q^ ^^q + h.c. (A.3) 

Note that a complex scalar in the A/" = 4 supermultiplet, as well as Mq, appear in the 
scalar terms after solving for the auxiliary field constraints within the full coupled theory. Note 
that both these operators have conformal dimension^^ A = 3, which matches the standard 
prescription for the powers of p appearing in the asymptotic behaviour of the fields - see 
below. 

Let us comment on the derivation of O^. The mass term for the hypermultiplet fields 
originates from the following superpotential term 

iV2 J (fe (q $7,7(3 -/i.e.) (A.4) 

^''This dimension applies in the UV where the effects of quark mass are neghgible and the theory becomes 
conformal. 
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where Q and Q are chiral superfields containing {q^il)) and (q,tp), respectively. These hy- 
permultiplet fields appear as ground states of the 3-7 and 7-3 strings while, as the subscript 
indicates, the superfield describes a particular set of massless modes in the 7-7 string sec- 
tor. In particular, the lowest component of is a complex scalar describing the transverse 
fluctuations of the D7-brane position, i.e., 

for the orientation in eq. (3.1). After the decoupling limit, this is no longer a dynamical 
field in the gauge theory but its expectation value sets the mass of the hypermultiplets, 
i.e., ($7,7) = Mq/\/2. One sees from eq. (A. 5) that the geometric angle (f) appearing in 
our construction of the D7-brane embeddings in section 3 corresponds precisely to the phase 
of the complex field $7,7. Hence in deriving O^, we consider a phase rotation with the 
given expectation value for $7^ in the Lagrangian (A.4). The resulting variation yields 
5C = 6(1) Mq and so we have divided by the factor of Mq in defining the operator given in 
eq. (A.3). 

Recall that the full flavour symmetry is U{Nf), which of course matches the worldvolume 
gauge symmetry of the D7-branes. The current is the conserved current corresponding to 
the diagonal L'"(l)q of this global symmetry, i.e., J* is the quark charge density. Our discussion 
can easily extended to the SU{Nf) symmetry by considering the nonAbelian gauge fields on 
the D7-brane. The corresponding fiavour currents would be 

^ = iT'% ' [r&^i^j + V^V^V;] - ^ [q^'V^qj {V^qfqj) (A.6) 

-i{qH'D^q)]-V^^rq})] , 

where T" are (Hermitian) generators of SU{Nf) and we have restored explicit flavour indices 
on the fields. With the operators and O^, we have also focussed on the SU{N() neutral 
terms but it would be straightforward to extend these to a more general discussion. For 
example, in general transforms in the adjoint of U (Nf) and so one can easily chose a 
more elaborate mass matrix rather than one proportional to the identity as above. 
Recall the equation of motion (3.7) for x implies that asymptotically 

X = - + ^ + --. (A.7) 

The dimensionless constants m and c are related to the quark mass and condensate as [10] : 

Mq = iVArm, (A.8) 

{OJi = ~V\N,N:r^c. (A.9) 

As presented here, this dictionary was established for a constant coefficient {i.e., the mass) and 
uniform expectation value of the operator O^. However, the same relationships also apply 
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in considering correlators where Mq is shifted with a general space- and time-dependent 
coefHcient or source n{x). In section 4.2, we express the relevant correlators in terms of 
variations S9{x,p) (or its Fourier transform along the x^^ directions). Given that x = cos^, 
we have SO ~ —Sx asymptotically where x approaches zero. Hence, as confirmed from 
eq. (4.44), 50 has the following the asymptotic behaviour: 

Se{x, p) = ^ ' + + ■■■ , (A.IO) 
P P"^ 

The source term is related to the first coefficient above as in eq. (A. 8) (up to a sign): ^{x) = 
— ^y/XTdOo{x). Similarly the induced expectation value and @{x) are related as in eq. (A. 9), 
again up to an overall sign. 

Ref. [38] investigated the present holographic theory at finite quark density and so estab- 
lished a similar dictionary for the (time component of the) worldvolume gauge field. With 
the asymptotic behaviour 

At = fi- - + ■■■, (A.ll) 

P 

it was found that fj, is precisely the chemical potential or the coefficient of the charge density 
operator and the quark density was given by 

ng = {Jl) = ^N,N,T^a. (A.12) 

Again these relations were originally established for a constant chemical potential and uniform 
quark density but they still apply for more general configurations. In particular for the 
correlators of section 4.1, we consider more general gauge field configurations with asymptotic 
behaviour 

A^ = S^(x) + ^ + ---. (A.13) 

In this case, T^i^ corresponds to the (space- and time-dependent) coefficient of the current 
operator and the induced expectation values are given by 

{J^i^)) = ^^f^c '7^{x) . (A.14) 

Note that here the index on a'^ is raised with jy'*^, the inverse metric in the CFT.-^^ 

Now we would like to turn to the holographic dictionary for the dcf) modes. Prom the 
equation of motion (4.43) and the asymptotic behaviour of % given in eq. (A. 7), we can 
determine that Scf) has the following asymptotic behaviour 

6(P{x,p) =6Mx) + ^^ + --- ■ (A.15) 
P^ 



^^Note that while the calculation of the spectral function in section 4.1 is presented in terms of the gauge- 
invariant field strengths this was simply choice of convenience and the final spectral function corresponds 
to that for current in eq. (A.l). 
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From eq. (A. 5), we saw that the geometric angle appearing in the D7-brane embeddings 
corresponds precisely to the phase of the complex field $7^. Hence S(I)q corresponds to 
a fluctuation in the phase of the hypermultiplet mass term. As usual, the dimensionless 
constant $ is proportional to the induced expectation value of but we would still need 
establish the precise constant of proportionality. 

To establish the latter relationship, it is natural to frame the discussion in terms of the 
thermal partition function - see [38], for example. -"^^ The source potential (50o enters the 
partition function as 



exp [-(3F{P, S(f)o)] = J2 exp j d^x {H - 6(j)o O^) 



(A.16) 



where a sum over all states is denoted on the right hand side and /3 denotes the inverse 
temperature. Of course, F{l3,6(f)Q) and H are the free energy and Hamiltonian densities, 
respectively. To begin, we note that, as can be seen from eq. (A.16), 



5F 



(A.17) 



To compare to the string description, we turn to the semiclassical analysis of the Euclidean 
supergravity path integral - see [10, 38], for example. Here the on-shell action gives the 
leading contribution to the free energy, i.e., 1^ = l3F. Hence to compare to eq. (A.17), we 

need to evaluate the change of the on-shell D7-brane action induced by a variation Scbo. The 
background solution for this field is simply cp = and so to linear order the action is invariant. 
Hence we can focus on the appropriate Euclidean version of the quadratic action (4.45) and 
the variation yields a boundary term, as in eq. (4.47), 



SF = ^N.T^rro^ J 



pV/(i-x')V 



p—*oo 



(A.18) 



where the last expression uses the asymptotic behaviour of both x ^4* from eqs. (A. 7) 
and (A. 15), respectively. Comparing eqs. (A.17) and (A.18), we find 



vr^ 1 



2M, 



■NtT^.TQ^m' ^{x) = -N,N, M^T' $(x) , 



(A.19) 



which completes the holographic dictionary for the S(f) modes. 

Above we have identified the i = modes of the worldvolume fields with operators in the 
dual field theory. Of course, a similar dictionary relates the higher-^ modes to dimension i + 3 
operators in the gauge theory. Qualitatively, we may say that the latter are constructed 



In principle then this discussion only concerns space- but not time-dependent sources. However, the 
following results apply for the general case including time dependence. 
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with products of i adjoint scalar fields 'sandwiched' between two fundamental fields in the 
above operators - see e.g., [11, 56]. For example, the expression (A.4) for would be 
generalized to 

j d^9Q^r,7 ($3,3)^ Q (A.20) 

where (<&3,3)^ represents a traceless combination of scalar superfields in the adjoint hyper- 
multiplet.^'' As is often the case in the AdS/CFT correspondence and its generalisations, 
the precise matching between normalisations in the field theory and string theory is difficult, 
mainly due to the fact that one would need the full D3/D7 brane action before taking the 
decoupling limit to determine the couplings of the bulk fields to field theory operators. 

However, to study the spectral functions overall numerical constants need not concern 
us and we simply choose a convenient normalisation which is consistent with holography. 
However, we must ensure that the spectral function has the proper scaling dimension (2.14). 
The three equations of motion, (4.19) for the vector, (4.53) for the pseudoscalar, and (4.59) 
for the scalar imply the asymptotic behaviour 

^t = Aep^ + Bip-^-\ (A.21) 

for some constants A^, where '^i = E^., S(l)£, p56i. In order to obtain the correct scaling 
of the spectral function, wc change to the standard dimensionful coordinates z = I? jr = 
\/2/7rTp used in the usual AdS/CFT prescriptions for computing correlation functions — see, 
e.g., [4]. In these coordinates, (A.21) becomes 

= Ae z'^ + Be _ ^^22) 

Hence, for z ^ (the boundary), wc expect the leading behaviour ^' ~ z~^. Taking a cutoff 
at small z = e, this implies that we should take 

^e{p,k)-^Ak)e 2^/2 ^^^p^y (A-23) 

where represents E£_k for the vector, Vg^k for the pseudoscalar, and pTZg^k for the scalar. 
As we will see in the next section, eq. (A. 23) ensures that the spectral functions have the 
correct asymptotic scaling (2.14) appropriate for an operator of dimension A = ^ + 3. 



B. Spectral function for scalar meson at T = 

In the low temperature phase of the D3/D7 brane theory (in the large- A^^ limit), we expect 
the scalar spectral function to be a series of (5-functions positioned at mass eigenvalues of 
the corresponding mesons. In this appendix, we explicitly demonstrate this expectation by 

^^Hcncc in this expression (A.20), the fundamental fields have an implicit sum of over the global U{Nf) 
indices with $7,7 and the gauge U{Nc) indices with ($3,3)*. 
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calculating the scalar spectral function at zero temperature. The mesons and their spectrum 
in this theory at T = were studied in detail in ref. [11]. 

The background geometry corresponding to the field theory at T = is that of (2.15) with 
ro = so that /(r) = 1, i.e., the background no longer contains a black hole. We embed A^f D7- 
branes in this geometry as described by the array (3.1), at a distance mo = 27r£^Mq from the 
D3-brancs. The resulting configuration is super symmetric with eight supercharges preserved. 
Using spherical polar coordinates in the 4567-space with radial coordinate g = moQ, the 
induced metric on the D7-branes is 

^(1 + g') i-dt' + dx2) + ^ {dg' + g'dnl) . (B.l) 



ds^ 



Following the same procedure as described in section 4.2, we consider small scalar fluctu- 
ations 5R of the D7-branes about this fiducial embedding. That is, taking polar coordinates 
(j) and R = moR in the 89-directions, we consider (f) = and R = mo(l + SR). Expanding the 
DEI action to quadratic order in 6R, we find 



The corresponding equation of motion is 

da 



+ Q- 



^g^'dbSR 



^g^'^daSRdbSR. 



0, 



(B.2) 



(B.3) 



where g is the induced metric (B.l). 

We take the fluctuations to be constant on the internal S^. Then integrating over the 
and integrating by parts, the action (B.2) becomes 



S = -TT^mt Tnr Nc J d'^x [g^5R dg5R] 



»oo 



(B.4) 



Expanding the fluctuation in terms of Fourier modes in the worldvolume directions, we take 
5R{k, g) = SRo{k)Tlk{g) /TZkiQoo) as the solution of eq. (B.3) regular at ^ = and normalized 
to SRo{k) at g = goo- The correlation function is given by 



5— >oo 



nk{g) 



(B.5) 



£1— »0O 



Changing the radial coordinate from g to z = g^/{l + g'^) and setting fc^ = (u;, 0,0,0), the 
equation of motion (B.3) for a fluctuation with vanishing spatial momentum can be written 
as 



4z(l 



-7^fc(^) = 0, 



(B.6) 



with UJ = CO L /rriQ. The solution regular at z = is given by a hypergeometric function 



Tlkiz) = 2F1 Q + ^Vl + a)^ I - \Vl + 0^; 2; z^ . 



(B.7) 
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In [11], the meson spectrum was determined by adjusting u) such the this radial wavefunction 
vanishes asymptotically as z ^ 1. Hence the correlator (B.5) has a series of poles positioned 
precisely at the mass eigenvalues of the meson spectrum. The imaginary part of the correlator 
will then come from deforming the integration contour for lo into the complex plane, in the 
usual way. The difference between various types of correlators {e.g., retarded, advanced, 
Feynman) then comes from the precise choice of this contour. 

With the solution in eq. (B.7), the correlator (B.5) becomes-^^ 



^ 2 7r2m^ro,iVf _2 

^ — TTo ^ 



, /I 1 A ^\ vr nVTTo^ 



(B.8) 



where we dropped contact terms. As indicated above, the correlator (B.8) has poles at 
o;^ = 4n(n + l)mQ/L'^, n = 1,2, corresponding to the meson spectrum found in [11]. Using 
the expansion 

OD 

IT vra 



71 /I u v-^ 

— tan — = > 

An o 



and Sokhotsky's formula 



we find 



4a 2 ^ (2n + 1)2 - 

n=0 ^ ' 



lim — - — = ^iTrS(x) +V {- 

e^o x±ie \x 



oo 



^ = -2Im G = ^ ^ yrr^J {u' - ul) . (B.9) 

TT ^ — ' 

n=l 

Since 

the spectral function for a; > can be expressed as 

AT AT °° I 1~ 

^ = ^Y.^U l + ^SicO- Un) . (B.ll) 

ri = l V " 

This expression confirms our expectations for the spectral function in the low temperature 
phase of the theory, illustrated in fig. 2a: the spectral function is a sum of (5-functions posi- 
tioned at the meson mass eigenvalues. Of course, the same is true in the vector and pseu- 
doscalar channels (as well as for ^ > in any of the channels). Here we have explicitly 
calculated the spectral function only for T = because we can analytically solve for the 
radial profile of the fluctuations. However, for the general case with < T < T{^^^, the calcu- 
lation is similar and the spectral functions again would take the form of a sum of (5-functions 
positioned at the meson masses. 

For large values of u (at which large n is relevant), the above result (B.ll) can be written 

as 

^ ~ AcOn S {U - Un) ■ (B.12) 



'in order to obtain the appropriate normalization we must multiply G by 1/Mq - see section 4.2. 
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where Aa)„ = (Dn+i — — 2 is the spacing between delta functions for large n. Note that 
the factor N^N^uP' /Air is precisely the high frequency asymptotics found for the £ = spectral 
functions in appendix C - see eq. (C.7). With the appropriate choice of coordinates, the 
radial equation in appendix C takes the same form as here, however, the boundary conditions 
chosen there for small radii assume a black hole embedding and naturally produce a complex 
result for the radial profile. This can be contrasted with the (naively) real result (B.7) for 
the Minkowski embeddings here. However, eq. (B.12) shows that if we examine the low 
temperature spectral function (B.ll) with low resolution on the scale of the spacing Aa)„, 
the (5-functions in the spectral function are averaged out and we reproduce the same high 
frequency tail as for the black hole embeddings. 



C. Spectral function high frequency asymptotics 

In this section we find expressions for the vector, pseudoscalar and scalar spectral functions 

in the high frequency limit, i.e., uo much larger than all scales: u) 2> T,M^. Note that in a 
sense this limit is equivalent to taking the limit of zero temperature and quark mass. 

The spectral functions, '^^^{u}), 'H^^(ct;), Tii^Lv), are collectively denoted by here 
and are defined as 



2£ 



NtN,T^^+^ Im 



(C.l) 



J p— >oo 



where = pTZi^ki 'Pe,k, E^^k for the scalar, pseudoscalar, and vector fluctuations, respectively. 
The desired limit is achieved by considering both p ^ oo and to^ — ^ +cxd. In this limit, 
the embedding function is x ~ m/p, and the three equations of motion, (4.26) for Ek, (4.59) 
for TZk, and (4.53) for Vk, reduce to 



3 
— ( 

P 



2) 8(«)2 - q 



2M 



(C.2) 



Changing variables first from p to the 'standard' radial coordinate r using (3.2) (which 
becomes r = pro/V2 asymptotically) and then to z = L? /r, we obtain 



$^'(z) - -$;(z) 



+ fe2U^(z) = 0, 



(C.3) 



where k = —u + q . For timelike momenta, the solution satisfying the incoming wave 
boundary condition at z = oo (the horizon in the zero temperature limit) can be written in 
terms of the Hankel function of the first kind. 



^ = zHf\{\k\z), 



(C.4) 



assuming u; > [15]. The spectral function then becomes 



27r2 



Im 



(C.5) 
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For £ = 0, this is 



47r 



which for vanishing three-momentum q = reduces to 

which coincides with the high frequency asymptotics which we found for all three £ = 
spectral functions in section 4. Of course, this asymptotic behaviour precisely matches that 
of the analytic solution (4.33) with £ = found for the vector modes in the massless quark 
limit. 

For £ = 1 and vanishing three-momentum q, the spectral function is 

(C-8) 

Again, this matches the asymptotic behaviour of the analytic vector solution (4.33) for Mq = 
with £ = 1. In fact, given that these asymptotics are independent of Mq and that we have 
found a common expression for all three channels, we can use eq. (4.33) to write the general 
asymptotic behaviour with g = as 

Hence we have produced the expected asymptotic behaviour (2.14) for an operator of scaling 
dimension A = £+3. Of course, the behaviour for general q is in principle given by combining 
the expressions (C.4) and (C.5). 



D. Effective potentials and quasinormal modes 

In this section, we rewrite the equations of motion for the pseudoscalar and scalar fluctuations 
of the D7-brane in the form of the Schroedinger equation. This can be considered a first step 
towards calculating the spectrum of quasinormal modes for these fields - see, e.g., [36]. The 
effective potential in each of these effective Schroedinger problems allows us to infer certain 
aspects of the quasinormal spectra. In particular, we argue that tachyonic modes appear in 
the scalar spectrum sufficiently close to the critical solution. The same analysis for the vector 
modes gives results that are essentially identical to those found for the pseudoscalar. 

D.l Pseudoscalar 

Considering fluctuations of the pseudoscalar, we take the general ansatz 

6cP^e"'^V{p)ye{S^), (D.l) 
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where yeiS"^) are spherical harmonics on the and = (—io,q,0,0). Then the pseu- 
doscalar's wave equation (4.43) can be written as 



dp [Hi dpP] + [q2 H2 + L^Hs]r = vo^r , 



with: 



(D.2) 



(D.3) 



/2 



Then, with the substitution V = hip, eq. (D.2) becomes 



L2 = £(^ + 2) 



q2 H2 + L^H3- Ho 



h Hih 



tp = xx> tp 



(D.4) 



where, as usual, the dot denotes a derivative with respect to p. The first term above can be 
rewritten as 

-Ho4> = - v^dp (yWodpip) + -H01P = -dl. iP + -Ho^ (D.5) 

where 

dp 



R* = 



(D.6) 



In terms of this new radial coordinate, the second derivative term takes the simple form found 
in a one-dimensional Schroedinger equation. Combining eqs. (D.4) and (D.5), all of the terms 
involving ip are eliminated if we choose h as 



1/4 



H 



1/2 



Hence the radial equation reduces to 



-dl,iP + Veffi; = SiP 
where the effective energy and potential are given by 



Veff=Cf''H2 + L^H3-Ho 



h Hi h 



(D.7) 



(D.^ 



(D.9) 



Let us comment on the new radial coordinate. In some sense, this coordinate is like 
the 'tortoise' radial coordinate introduced in the analysis of physics in the Schwarzschild 
geometry [57]. Approaching the event horizon, i.e., as /> — 1, i?o — (p — 1)^ and so R* oc 
— log(p — 1) +00. For large p, Hq ~ and so R* ~ 1/p 0. Note that given the 
definition in eq. (D.3), Hq is positive everywhere on the range p G (l,oo). Hence we are 
assured that R* is a monotonic function of p. 
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Although R* is the appropriate coordinate to analyse the effective Schroedinger equation 
(D.8), we can still gain some intuition for the problem by plotting the effective potential (D.9) 
as a function of p for various different D7-brane embeddings. A summary of results (with 
q,i = 0) is given in fig. 15. Note that for all embeddings, the effective potential exhibits 
a large barrier in the asymptotic region, as expected for an asymptotically AdS geometry. 
Further, in all cases, the effective potential vanishes at the horizon. For small xoi the poten- 
tial is monotonically increasing with p. For larger xoi & small potential barrier develops at 
intermediate values of the radius. Intuitively, the latter might give rise to metastable states 
in the effective Schroedinger problem. 




Figure 15: The effective potential (for 5(j) fluctuations) versus p for various xo with £ ~ and q — Q. 

For i > nonzero and q ^ the results are qualitatively similar to those depicted in 
fig. 15. The effective potential vanishes at the horizon and grows as p ^ oo. The potential 
barrier grows most quickly for £ ^ 0, due to the term proportional to in (D.9) which is 
roughly i{i + 2)p'^ /8 for large p. For values of xo near 1, a small potential barrier develops for 
intermediate values of the radius. Note that while a small potential barrier is already evident 
for Xo = 0.9 (with £ = and g = 0) in fig. 15, the potential barrier only develops for larger 
values of xo; Xo > 0.99 for £ > and/or q 0. 

The quasinormal modes of the pseudoscalar can be found by solving the Schroedinger 
problem constructed above, with the appropriate boundary conditions. One of the boundary 
conditions is that the pseudoscalar should have only an incoming wave component at the 
horizon, i.e., p — > 1 or i2* — > oo. Given the ansatz (D.l), we are thus looking for solutions 
with^^ ^ oc exp(zt»i?*). For large p or small R* where the effective potential diverges, 
we demand that the wavefunction vanish. Generically, these boundary conditions lead to 
complex eigenvalues for the effective energy £, which is in accord with our expectation that 
the quasinormal frequencies have the form [58, 59] 

m = ±n-iT, (D.IO) 

^^Note that taking the hmit p ^ 1 in eq. (D.7) yields a simple constant for h and so we have V cc tp as we 
approach the horizon. 
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with ri, r > 0. Note that F > ensures that the quasinormal excitations decay in time, as 
can be seen from the ansatz (D.l). However, given that 

£' = ro2 = (i^2_p2^-p2if^r, (D.ll) 

some translation is required to use our intuition for the Schroedinger problem to infer general 
characteristics of the quasinormal spectrum. Note that as the sign of Im{S) is not fixed, 
we are implicitly admitting energy eigenvalues which would correspond to both decaying and 
growing wavefunctions in the effective Schroedinger problem. However, this Schroedinger 
problem is purely an auxiliary tool and so one should not ascribe any physical significance to 
this observation. 

Above, we observed that for small XO) the effective potential rises monotonically from zero 
as we move away from the horizon towards larger p. Hence we would infer that Re{£) > 
or > r. Further we should not expect any suppression of Im{E), i.e., Im{£) ~ Re{£), 
which means that we should still expect to be ft and F to be the same order of magnitude in 
this regime. This intuition would then suggest the absence of any interesting structure in the 
corresponding spectral functions, as the quasinormal frequencies should be far from the real 
axis. However, as noted above, a small potential barrier appears at intermediate values of R* 
(or rho) as xo approaches one. Intuitively, then one might expect to find long-lived states with 
Re{£) Im{£) and Re{£) ^ 14//(well), i.e., Re{£) would be roughly given by the height of 
the potential in this intermediate potential well. From eq. (D.ll), this requires Q F with Q 
finite and so would correspond to quasinormal frequencies approaching the real axis. Hence 
we would expect the formation of peaks in the spectral function in this regime, as discussed 
in section 2. Of course, this intuitive picture developed from the effective potential matches 
the behaviour of the spectral functions found in section 4.2.1. We emphasize, however, that 
this intuition only gives a very rough picture of the quasinormal spectrum and it would be 
interesting to develop more detailed picture with a full calculation. 

As mentioned above, the results for the effective potential for the vector fluctuations are 
essentially the same as for the pseudoscalar and hence the quasinormal spectrum should also 
be similar. We saw in section 4.1 that the behaviour of the vector spectral functions is very 
similar to those for the pseudoscalar. 



D.2 Scalar 

With the ansatz 69 ~ e**^^ 7?.(p) 3^^(S'^)), the scalar wave equation (4.44) can be written as 
Ho 



HidpTZ 



+ [Ki^ H2 + {i + ^){i-i)H^]n = w'^n 



(D.12) 



where i?o, H2 and are the same as defined in eq. (D.3), while we must redefine the 
following: 



^^-(r?f^7^'^^-(^+3)(.-i). 



(D.13) 
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Since eq. (D.12) has the same form as eq. (D.2), we can use precisely the same steps 
as above to cast this equation for the scalar fluctuations into the form of a Schroedinger 
equation. That is, taking IZ = hifj with h = Hq^^ / h\^'^ and defining the radial coordinate 
(D.6), eq. (D.12) takes the form of the Schroedinger equation (D.8) with the effective energy 
and potential given by (D.9). Examining the effective potential to gain some intuition for the 
physics, we find: Again for all embeddings, there is a large potential barrier in the asymptotic 
region and the effective potential vanishes at the horizon. For embeddings of the D7-branes 
with < xo < 0-7, fig- 16 shows that the effective potential is a monotonically increasing 
function of p. For I = Q and (? = 0, once xo ^ 0.7, the potential develops a negative well near 
the horizon. As xo increases towards 1, this well near p ~ 1 becomes deeper and wider. Note 
that Xo — 0.9621 and xo — 0.99973885 correspond to the first and second kinks, respectively, 
in a plot of the free energy versus temperature - see fig. 4. For any modes with £ > 0, one 
finds that there is never a region where the effective potential becomes negative, however, for 
Xo near 1 the potential develops a small barrier near the horizon. For nonvanishing spatial 
momentum 7^ 0), the effective potential exhibits a negative well near the horizon for values 
of Xo near 1. However, the well is neither as deep nor as wide as that for g = 0. 

6 . ^ ^^^^^^^^^^^^^ 

2 3 4 5 6 7 8 

P 




Figure 16: The effective potential for the scalar field [50) versus p for various xo with t = {) and 
A; = 0. 

Certainly, the most interesting feature of the effective potential for the scalar is the 
negative potential well which develops and grows as xo — ^ 1- One would expect that if this 
negative well grows large enough, it will be able to support a 'bound' state with <S < 0. 
Actually, since such a state would still see a finite potential barrier between the center of the 
well and the horizon, it would still be a long-lived state. Using a WKB approximation, we 
can estimate that a (zero-energy) bound state will appear for [60, 61] 



1 

n I vr 

2 



r dR*J-V,ff{R*) (D.14) 
JRn ^ 



Ro 
P° dp 



'1 V-nO * 

where n is a positive integer and the integration is over the values of p for which the potential 
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is negative. A plot of //vr + 1/2 is given in figure 17 (for i,q = 0). This quantity reaches 1 
for Xo — 0.98297 and 2 for xo — 0.99986, and so we expect that the first two bound states 
form at roughly these values of xo- Below we will argue that the appearance of these bound 
states can be associated with the dramatic spikes that were observed in the scalar spectral 
functions in section 4.2.2. Not then that our results here do not (quite) match the values of 
Xo corresponding to the first two kinks (on the black hole branch) of the free energy - recall 
that the latter correspond to xo = -9621 and .99973885, respectively. However, we expect 
this discrepancy is likely due to the approximations inherent in the WKB calculation. 

2.5 

2 .25 

2 

1 . 75 

n 1 

jT^y 1.5 
1.25 
1 

. 75 

0.9 0.92 0.94 0.96 0.98 1 

Xo 

Figure 17: Plot of I/tt + 1/2 versus xo for £,q = 0. 

As the effective potential is positive and monotonically increasing in the small xo regime, 
we expect the eigenfrequencies (D.IO) in quasinormal spectrum will again have 0, and T with 
the same order of magnitude. Of course, the regime with large xo is more interesting because 
of the appearance of bound states. These modes with Re{£) < are distinguished since > 
0^, as seen in eq. (D.ll). Further as we alluded to above, the corresponding wavefunctions 
are below a potential barrier as i?* — > oo and so must have the form ip ~ exp(— |r| R*) to 
avoid a divergence at the horizon. Given the boundary condition there, i.e., ip oc exp{im R*), 
this requires that F < for these modes. Further then, it follows that these exceptional modes 
grow rather than decay in time and so these bound states really represent an instability of 
the system. 

The above discussion indicates that these bound states appear when a quasinormal fre- 
quency crosses the real axis and so their appearance should be signalled by a pole appearing 
in the scalar spectral function calculated in section 4.2.2. Further, however, we argued that 
as the eigenfrequency crosses the real axis, it moves from a regime where > F^ for F > 
to Q2 < F^ for F < 0. Hence at the point that F = 0, we must also have Q = 0.'^^ Hence 
we see that the quasinormal frequencies must be cross the real axis by passing through the 
origin. This is, of course, precisely what was observed in section 4.2.2, where the poles in 

^"Essentially this we are just saying that Re{£) — just as the bound states form. 
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the spectral function appeared at precisely a; = 0. Further, the lack of much structure in 
the spectral function aside from these poles would indicate that the quasinormal frequencies 
approach the origin uniformly so that we never find eigenfreqencies with O S> F. We reiterate 
that this discussion only gives a schematic picture of the quasinormal spectrum and it would 
be interesting to develop a more detailed picture with a full calculation [36]. 

E. Diffusion constants for Dp/Dq systems 

This appendix extends the computation of the diffusion constant using the membrane paradigm 
[17] described in section 5.1 to that for the gauge theory dual to the supergravity configuration 
of a Dq-brane probe in the near-horizon black Dp-branc geometry. 

The background geometry (2.15) is generalised to the near-horizon black Dp-brane metric 
(in the string-frame) [3]: 

ds-" = H-l {-fde + d^l) rH^l_)j , (E.l) 

where now H{r) = {L/ry^'^ and /(r) = 1 — (ro/r)^^*'. The background also includes non- 
trivial dilaton and RR fields: e'^ = C'oi...^ = . The Hawking temperature 
associated with the horizon at r = ro is given by 



-P ro\^ 



AttL \L 



(E.2) 



According to the gauge/gravity correspondence, string theory on this background is dual to 
a supersymmetric (p -|- l)-dimensional gauge theory at temperature T. 

Consider placing a probe Dq-brane in the above geometry such that the probe has d 
spatial directions parallel and n + 1 transverse to the background Dp-branes, so that q = 
d+ n + 2 and such that it intersects the horizon at r = tq. In analogy to eq. (3.2), it is useful 
to introduce a new (dimensionless) radial coordinate p related to r via 



(rop) = r (^-f + 7-P - r^'f . (E.3) 

The horizon is now positioned at p = 1. Implicitly, we will assume in the following that the 
Dp/Dq system under consideration is T-dual to the D3/D7 one described by the array (3.1). 
This choice ensures that the brane configuration is supersymmetric at zero temperature and 
the probe brane embeddings should be stable in the finite temperature background (E.l) [10]. 
With the coordinate (E.3), the metric and the dilaton may be written as: 



2 1 /rop\(7-P)/2 
= 2 [^) 



P2 



-Ldt^+fdxi 

(p-3)/2 

/rop\ (7-p)(p-3)/4 



-I- h{p) [dp'^ + p^ [de'^ + sin^ Od^ln + cos^ ed^-j-p-n)] 
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where 



VP 



7-p 



Kp) = Vp 



7-p 



Hp) = rliL/rop) 



(7-p)/2 



(p-3)/(7-p) 



Describing the probe brane profile using xip) = cos6{p), the induced metric on the 
Dq-brane may be written as ds'^{g) = ds^{g) + Z{p)d^^, where 



1 /ro/)\(7-p)/2 

2 VT) 



ds\g) 

Z, \ 1 

Z{p) = h{p)p\l-x^ 



+ Kp) 



1 - X + P X 



■dp\ 



Using the DBI action and expanding the gauge fields to quadratic order, the relevant portion 
of the action for the gauge fields is 



I<i,F = -Tgiiri'tyn^ I dt d^x dp V^F^, g\ 



5, 



eff 



eff 



(E.4) 



We are now in a position to evaluate the diffusion constant using eq. (2.27) from [17]: 

-2 



D 



-g 



p=i 



~ , -gttgppgeff 

dp , 

1 v-f 



27rr 



^ P p {I- ^2)(n+l)/2 



(E.5) 

(E.6) 



where 

d — p (n — l)(p — 3) 



a = 



+ 



2(7 -p) 



(7-p)(p + n-3-d) 
(5 = n. 



_ (n- l)(p-3) d + A-p 
^~ 2(7-p) + 2 • 



One may check that foi p = 3 = n = d this result reduces to that for the D3/D7 case given 
in eq. (5.1). 

We have also evaluated DT numerically for the D4/D6 case {p = A, n = 2, d = 3) 
and the results are plotted in figure 18. The horizontal axis is labelled by the ratio of the 
temperature to the natural mass scale of the problem: 



M = 



AttL 



9eff{Mg) 



(E.7) 



Asymptotically, DT approaches 3/ An for large temperatures. As the temperature is reduced, 
DT decreases dramatically near the phase transition. The value at the phase transition is 
DT = .125 ~ .785/27r. If we continue following the black hole embeddings beyond the phase 
transition, DT continues to fall and it also becomes a multi-valued function of temperature, 
as was seen in fig. 13 for the D3/D7 system. Again, this simply reflects the fact that multiple 
embeddings can be found for a single temperature in the vicinity of the critical solution. 
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Figure 18: The diffusion constant D times the temperature T versus temperature T/M for a D6- 
brane probe in the black D4-brane geometry. For large T we have DT ~ 3/47r. 

Note that the asymptotic value for DT does not match that for the i?-charge diffusion 
constant calculated for a near-extremal D4-brane result: S/Svr [17]. However, there is no 
reason that these two quantities should be equal since the D6-brane does not fill the entire 
D4-brane throat, i.e., one of the D4-brane worldvolume directions is transverse to the D6 
probe. 

The D4/D6-brane system considered above is the basis for the construction of one holo- 
graphic model which mimics QCD at large [50]. Another interesting holographic model of 
a QCD-like theory comes from introducing D8 and anti-D8 probe branes in a D4-brane back- 
ground [51]. This system displays an interesting phase transition related to chiral symmetry 
breaking [62].^^ One can again calculate the diffusion constant for the quark charge in the 
high temperature phase along the lines described above. In this case, the D8-brane wraps the 
entire of the D4 background but otherwise fills the same directions as the D6-branes above. 
After the phase transition the embeddings are much simpler since there is no non-trivial ra- 
dial profile to be considered and, further, the embedding is temperature independent. The 
diffusion constant may be determined from eq. (E.6) setting x = = XO) p = 4 = n and 
d = 3. The result for the calculation is DT = 1/271. 
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